Precise Measurements of Direct CP Violation, CPT Symmetry, and Other 
Parameters in the Neutral Kaon System 
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We present precise tests of CP and CPT symmetry based on the full dataset of K — ¥ nn decays 
collected by the KTeV experiment at Fermi National Accelerator Laboratory during 1996, 1997, and 
1999. This dataset contains 16 million K — > 7r°7r° and 69 million K — > 7r + 7r~ decays. We measure 
the direct CP violation parameter Re(e'/e) = (19.2 ±2.1) x 10~ 4 . We find the Kl-Ks mass difference 
Am = (5270 ± 12)xl0 6 hs' 1 and the K s lifetime r S = (89.62 ± 0.05) xl0~ 12 s. We also measure 
several parameters that test CPT invariance. We find the difference between the phase of the indirect 
CP violation parameter, e, and the superweak phase, 4> t — cf>sw = (0.40 ± 0.56)°. We measure the 
difference of the relative phases between the CP violating and CP conserving decay amplitudes for 
K — s- 7r + 7r _ (4>+-) and for K —>■ ty°tt° (<poo), = (0.30 ± 0.35)°. From these phase measurements, 
we place a limit on the mass difference between K° and K°, AM < 4.8 x 10 -19 GeV/c 2 at 95% C.L. 
These results are consistent with those of other experiments, our own earlier measurements, and 
CPT symmetry. 

PACS numbers: 11.30.Er, 13.25.Es, 14.40.Df 



INTRODUCTION 



the 1964 discovery of CP violation in 
7T~ decay significant experimental effort has 



Since 

been devoted to understanding the mechanism of CP vi- 
olation. Early experiments showed that the observed ef- 
fect was due mostly to a small asymmetry between the 
K° — >• K° and K° — > K a transition rates, which is 
referred to as indirect CP violation. Decades of addi- 
tional effort were required to demonstrate the existence 
of direct CP violation in a decay amplitude. This paper 
reports the final measurement of direct CP violation by 
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the KTeV Experiment (E832) at Fermilab. 

Direct CP violation can be detected by comparing the 
level of CP violation for different decay modes. The pa- 
rameters e and e' are related to the ratio of CP violating 
to CP conserving decay amplitudes for K — » 7r + 7r~ and 



VOQ 



A(K L -> TT+TT") 
~ A(K S 7T + 7T~) 

A (K L -» A") 
= A(K S -> 7r°7r°) 



2e', 



(1) 



where e is a measure of indirect CP violation, which is 
common to all decay modes. The relation among the 

complex parameters rj^ , 7700, e, and e' is illustrated in 

Fig.ffl 

If CPT symmetry holds, the phase of e is equal to the 
"superweak" phase: 



<psw = tan -1 (2Am/Ar) . 



(2) 
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FIG. 1: Diagram of CP violating kaon parameters. For this 
illustration, the e parameter has the central value measured by 
KTeV and the value of e is scaled by a factor of 50. Although 
they appear distinct in this diagram, note that 4>+- an d 4>oo 
are consistent with each other within experimental errors. 



where Am = m^ — m$ is the K^-Kg mass difference and 
Ar = Y s — Ti, is the difference in the decay widths. 

The quantity e' is a measure of direct CP violation, 
which contributes differently to the 7r + 7r~ and ir°ir° de- 
cay modes, and is proportional to the difference be- 
tween the decay amplitudes for K° —> tt + tt~ (tt°tt°) and 
K° tt + tt~ (7r°7r°). Measurements of tttt phase shifts 
show that, in the absence of CPT violation, the phase of 
e' is approximately equal to that of e. Therefore, Re(e'/e) 
is a measure of direct CP violation and Im(e'/e) is a mea- 
sure of CPT violation. 

Experimentally, Re(e'/e) is determined from the double 
ratio of the two pion decay rates of Kl and Kg'- 

T(K L -> 7T+7T-) /T{K S ->• 7T+7T-) 



T(K L -> ttV) /T(K s -> 7r°7r°) 

2 

^tz. ttl + 6R e (e'/e). 
Voo 



(3) 



For small |e'/e|, Im(e'/e) is related to the phases of 774 
and 7700 by 



4>+- s» (f> e + Im(e'/e), 

(froo ~ (j>€ - 2Im(e'/e), 

A(j> = 0oo - ~ -3Im(e'/e) 



(4) 



The Standard Model accommodates both direct and 
indirect CP violation . Most recent Standard Model 
predictions for i?e(e'/e) are less than 30 x 10~ 4 
however, there are large hadronic uncertainties in these 
calculations. Experimental results have established that 



Re(e'/e) is non-zero [16H20|. The previous result from 
KTeV, which was based on about half of the KTeV 
dataset, is Re{e'/e) = (20.7 ±2.8) x 10" 4 [20[. This result 
was published in 2003 and will be referred to in this text 
as "KTeV03." The result based on all data from NA48 
at CERN is Re(e'/e) = (14.7 ± 2.2) x 10" 4 [ii]. 

This paper reports the final measurement of Re(e'/e) 
by KTeV. The measurement is based on 85 million re- 
constructed K — > tttt decays collected in 1996, 1997, and 
1999. This full sample is two times larger than, and 
contains, the sample on which the KTeV03 results are 
based. We also present measurements of the kaon pa- 
rameters Am and rg, and tests of CPT symmetry based 
on measurements of A<fi and 4> e — 4>sw ■ Using our phase 
measurements, we place a limit on the mass difference 
between K° and K°. 

For this analysis we have made significant improve- 
ments to the data analysis and the Monte Carlo sim- 
ulation. The full dataset, including the data used in 
KTeV03, has been reanalyzed using the improved re- 
construction and simulation. These results supersede 
the previously published KTeV03 results (20|, which were 
based on data from 1996 and 1997. 

This paper describes the KTeV experiment in Sec. UH 
the analysis technique in Sec. IIII1 and the extraction 
of physics results in Sec. IIVI We emphasize changes 
and improvements since the KTeV03 publication. We 
will refer to [13] for some details that have not changed 
since KTeV03. Section [V] presents the final KTeV re- 
sults, including correlations between the parameters and 
crosschecks of the results. Section IVII is a summary and 
discussion of the results. Appendix [Al contains a discus- 
sion of the dependence of our measurements on details of 
kaon regeneration. 



II. MEASUREMENT TECHNIQUE AND 
APPARATUS 



A. Overview 

The measurement of Re(e'/e) requires a source of Kl 
and Ks decays, and a detector to reconstruct the charged 
(7r + 7r~) and neutral (7r 7r°) final states. The strategy 
of the KTeV experiment is to produce two identical Kl 
beams, and then to pass one of the beams through a "re- 
generator" that is about two hadronic interaction lengths 
long. The beam that passes through the regenerator is 
called the regenerator beam, and the other beam is called 
the vacuum beam. The regenerator creates a coherent 
\K£)+p \Ks) state, where p, the regeneration amplitude, 
is a physical property of the regenerator. The regener- 
ator is designed such that most of the K tttt decays 
downstream of the regenerator are from the Ks compo- 
nent. The charged spectrometer is the primary detector 
for reconstructing K — > tt + tt~ decays and the pure Ce- 
sium Iodide (Csl) calorimeter is used to reconstruct the 
four photons from K — > 7r°7r° decays. A Monte Carlo 
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simulation is used to correct for the average acceptance 
difference between K — > tttt decays in the two beams, 
which results from the very different Kl and Kg life- 
times. The decay- vertex distributions provide a critical 
check of the simulation. The measured quantities are 
the vacuum-to-regenerator "single ratios" for K — > tt + tt~ 
and K — > 7r°7r° decay rates. These single ratios are pro- 
portional to |t7_| / p\ 2 and |?7oo/p| 2 , respectively, and the 

ratio of these two quantities gives Re(e'/e) via Eq. [3] 



B. KTeV Experiment 

The KTeV kaon beams are produced by a beamline 
of magnets, absorbers, and collimators that act on the 
products of a proton beam incident on a fixed target. 
The 800 GeV/c proton beam, provided by the Fermilab 
Tevatron, has a 53 MHz RF structure so that the pro- 
tons arrive in ~1 ns wide "buckets" at 19 ns intervals. 
This beam is incident on a beryllium oxide (BeO) target 
that is about one proton interaction length long. Im- 
mediately downstream of the target, the beam consists 
of protons, muons, and other charged particles, neutral 
kaons, neutrons, photons, and hyperons. This beam is 
collimatcd into two beams and the non-kaon component 
is reduced by magnets and absorbers in a 100 meter long 
beamline. At the start of the fiducial decay region, 120 m 
downstream of the target, the average kaon momentum 
is about 70 GeV/c. The neutron-to-kaon ratio is 1.3 in 
the vacuum beam and 0.8 in the regenerator beam. The 
KTeV beams and the beamline elements that produce 
them are described in detail in [20l |. 

KTeV reconstructs kaon decays that occur in an evac- 
uated decay region 90-160 m downstream of the target. 
Figure[2]is a schematic of the detector. In the KTeV coor- 
dinate system, the positive a;-axis points to the left if the 
observer is facing downstream, the positive y-axis points 
up, and the positive z-axis points downstream from the 
target. At the upstream end of the decay region, the re- 
generator alternates between the two beams to minimize 
acceptance differences between decays in the vacuum and 
regenerator beams. The charged spectrometer and Csl 
calorimeter are located downstream of the vacuum win- 
dow at the end of the decay region. The decay region and 
primary detectors are surrounded by a system of photon 
veto detectors to detect particles with trajectories that 
miss the Csl calorimeter. The major detector elements 
are described in more detail in the following paragraphs. 

The regenerator consists mainly of 84 10 x 10 x 2 cm 3 
scintillator modules as seen in Fig. [3^. Its primary pur- 
pose is to provide K$ regeneration, but it is also used 
as part of the trigger and veto systems. Each module is 
viewed by two photomultiplier tubes (PMTs), one from 
above and one from below. The downstream end of the 
regenerator has a lead-scintillator sandwich called the 
"regenerator Pb module" (Fig.O)), which is also viewed 
by two PMTs. This last module of the regenerator is 
used to define a sharp upstream edge for the kaon decay 



Regenerator 



Trigger Muon 
Drift Hodoscope Veto 
Chambers ^ CA 

Csl 



Beams 



. Regenerator Beam 
■ Vacuum Beam 



Mask Anti 



25 cm I 



Photon Vetoes 



_ Photon 
Analysis Vetoes 
Magnet Steel 



120 140 160 180 

Z = Distance from kaon production target (meters) 



FIG. 2: Schematic of the KTeV detector. Note that the ver- 
tical and horizontal scales are different. 



region in the regenerator beam. 

The charged spectrometer consists of four drift cham- 
bers (DCs) and a dipole magnet. Each drift chamber 
measures charged-particle positions in both the x and 
y views. A chamber consists of two planes of horizontal 
wires to measure y hit coordinates, and two planes of ver- 
tical wires to measure x hit coordinates; the two s-planes 
and the two y-planes are offset to resolve position ambi- 
guities. The DC planes have a hexagonal cell geometry 
formed by six field-shaping wires surrounding one sense 
wire (Fig. @}. There are a total of 1972 sense wires in the 
four drift chambers. The cells are 6.35 mm wide, and the 
drift velocity is about 50/xm/ns. The analyzing magnet 
imparts a kick of 412 MeV/c in the horizonal plane. The 
well-known kaon mass is used to set the momentum scale 
with 10 precision. 

The Csl calorimeter consists of 3100 pure Csl crys- 
tals viewed by photomultiplier tubes. The layout of the 
1.9 x 1.9 m 2 calorimeter is shown in Fig. [5j There are 
2232 2.5 x 2.5 cm 2 crystals in the central region, and 
868 5 x 5 cm 2 crystals surrounding the smaller crys- 
tals. The crystals are all 50 cm (27 radiation lengths) 
long. Each crystal is wrapped in 12 pm, partially- 
blackened, aluminized mylar in a manner designed to 
make the longitudinal response of each crystal as uniform 
as possible. The calorimeter is read out by custom dig- 
itizing electronics (DPMTs) placed directly behind the 
PMTs|2l|. Momentum- analyzed electrons and positrons 
from Kl — > ir e^v {Kes) decays are used to calibrate 
the Csl energy scale to 0.02%. 

An extensive veto system is used to reject events com- 
ing from interactions in the regenerator, and to reduce 
background from kaon decays into non-7T7r final states 
such as Kl — > 7r ± /i =F j/ and Kl — > 7r 7r°7r°. The veto sys- 
tem consists of a number of lead-scintillator detectors in 
and around the primary detectors. 

KTeV uses a three-level trigger to select events. Level 1 
uses fast signals from the detector and introduces no 
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FIG. 3: Diagram of the regenerator, (a) Layout of the 85 re- 
generator modules, including the lead-scintillator module, (b) 
Zoomed diagram of the lead-scintillator regenerator module. 
The PMTs above and below are not shown. The thickness of 
each lead (scintillator) piece is 5.6 (4.0) mm. The transverse 
dimension is 100 mm, and is not drawn with the same scale 
as the z-axis. The kaon beam enters from the left. The ar- 
rows indicate the location and ±lcr uncertainty of the effective 
upstream edges for reconstructed K — > and K — ¥ tt + -k 

decays for 1999 data. 



deadtime. Level 2 is based on more sophisticated pro- 
cessing from custom electronics and introduces a dead- 
time of 2-3 fj,s; when an event passes Level 2 the entire 
detector is read out with an average deadtime of 15 /is. 
Level 3 is a software filter; the processors have enough 
memory that no further deadtime is introduced. 

Individual triggers are defined to select K —> ir + n~ 
and K —> 7r°7r° decays; the trigger efficiencies are studied 
using decays collected in separate minimum-bias triggers. 
Additional triggers select decays such as Kl ~ > ir^e^v 



and Kl — > ir°ir°ir which are used for calibration and ac- 
ceptance studies. The "accidental" trigger uses a set of 
counters near the target to collect events based on pri- 
mary beam activity; these events are uncorrelated with 
detector signals that come from the beam particles and 
are used to model the effects of intensity-dependent ac- 
cidental activity. 

Several changes were made to the KTeV experiment to 
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FIG. 4: Diagram of drift chamber geometry showing six field 
wires (open circles) around each sense wire (solid dots). The 
solid lines illustrate the hexagonal cell geometry; they do not 
represent any physical detector element. The vertical dashed 
lines are separated by 6.35 mm and are used to define the 
track separation cut described in Sec. HUB 21 
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FIG. 5: Beamline view of the KTeV Csl calorimeter, showing 
the 868 larger outer crystals and the 2232 smaller inner crys- 
tals. Each beam hole size is 15 x 15 cm 2 and the two beam 
hole centers are separated by 0.3 m. The positive z direction 
is into the page. 
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improve data collection efficiency for the 1999 run. 

1. Neutral Beams. The proton extraction cycle of the 
Tevatron was improved from 20 second extractions 
(or "spills" ) every 60 seconds in 1996 and 1997 to 40 
second extractions every 80 seconds in 1999. The 
maximum available intensity was ^2 x 10 11 protons 
per second. In 1999, KTeV chose to take about half 
of the data at an average intensity of ~ 1.6 x 10 11 
protons/s and half at a lower average intensity of 
~ 1 X 10 11 protons/s as a systematic cross-check. 

2. Csl Calorimeter Electronics. During 1996 and 
1997 data taking, individual channels of the custom 
readout electronics for the Csl calorimeter failed 
occasionally. These failures account for half of the 
20% data-taking inefficiency during 1996 and 1997. 
They also affect the data quality and complicate 
the calibration of the calorimeter. All of the cus- 
tom electronics were re-fabricated and installed in 
the Csl calorimeter in preparation for the 1999 run. 
The re-fabrication of the chips was successful; no 
Csl calorimeter electronics had to be replaced dur- 
ing the 1999 run. 

3. Drift Chambers. The drift chambers required some 
repair due to radiation damage sustained during 
data taking in 1996 and 1997. About half of one 
drift chamber was restrung and a second chamber 
was cleaned. The drift chamber readout electronics 
were modified to allow the system to run at higher 
gain without causing the system to oscillate or trig- 
ger on noise. 

4. Helium Bags. Helium bags are placed between the 
drift chambers to minimize the matter seen by the 
neutral beams after leaving the vacuum decay re- 
gion and to reduce multiple scattering of charged 
particles. In 1996 and 1997, one of the small he- 
lium bags was leaky and contained mostly air by 
the end of the 1997 run, so it was necessary in the 
analysis to correct for the increased multiple scat- 
tering resulting from this increased material in the 
detector. The bags were replaced for 1999. This 
reduction in material traversed by the beams was 
offset by a change in the buffer gas used in the drift 
chambers; the total ionization loss upstream of the 
Csl calorimeter was less than 5 MeV in each year. 
This energy loss occurs mostly in the scintillator 
hodoscope just upstream of the Csl calorimeter. 

5. Trigger. In 1999, the trigger was adjusted to select 
more Kt ~> 7r + 7r~7r decays for the measurement 
of kaon flux attenuation in the regenerator beam, 
called "regenerator transmission." The improve- 
ment of this measurement reduces several system- 
atic uncertainties associated with the fitting proce- 
dure as described in Sec. IIVB1 



C. Monte Carlo Simulation 

KTeV uses a Monte Carlo (MC) simulation to calculate 
the detector acceptance and to model background to the 
signal modes. The different Kl and K$ lifetimes lead to 
different z-vertex distributions in the vacuum and regen- 
erator beams. We determine the detector acceptance as 
a function of kaon decay z-vertex and energy, including 
the effects of geometry, detector response, and resolu- 
tion. To help verify the accuracy of the MC simulation, 
we collect and study decay modes with approximately ten 
times higher statistics than the K — > nir signal samples: 
K L -t n ± e T v, K L -> 7T 7r 7r°, and K L -t 7T + 7r~7r°. 

The Monte Carlo simulates K° j K° generation at the 
BeO target following the parameterization in 22] , prop- 
agates the coherent K° /K° state through the absorbers 
and collimators along the beamline to the decay point, 
simulates the decay including decays inside the regen- 
erator, traces the decay products through the detector, 
and simulates the detector response including the digi- 
tization of the detector signals and the trigger selection. 
The parameters of the detector geometry are based both 
on data and survey measurements. Many aspects of the 
tracing and detector response are based on samples of 
detector responses, called "libraries," that are generated 
with GEANT[23| simulations; the use of libraries keeps 
the MC relatively fast. The effects of accidental activity 
are included in the simulation by overlaying data events 
from the accidental trigger onto the simulated events. 
The Monte Carlo event format is identical to data and 
the events are reconstructed and analyzed in the same 
manner as data. More details of the simulation are avail- 
able in |20j . 

Many improvements have been made to the MC simu- 
lation since KTeV03[20j. We have improved the simula- 
tion to include finer details of electromagnetic showering 
in the Csl calorimeter and charged particle propagation 
through the detector. These changes are described in 
detail below. 

1. Shower library. For this analysis, the GEANT- 
based library used to simulate photon and electron 
showers in the Csl calorimeter has been improved 
to simulate the effects of incident particle angle. 
The library used for KTeV03 was binned in en- 
ergy and incident position. There were 325 posi- 
tion bins and six logarithmic energy bins (2 GeV, 
4 GeV, 8 GeV, 16 GeV, 32 GeV, and 64 GeV). 
The effect of angles was approximated by shifting 
the incident position based on the angle of inci- 
dence. The shower library has now been expanded 
to include nine angles in x and y (-35 mrad to 35 
mrad) for photons and 15 angles in x and y (-85 
mrad to 85 mrad) for electrons. Electron angles 
may be larger than photon angles because of the 
momentum kick imparted by the analyzing magnet. 
The position and energy binning is unchanged from 
KTeV03. Differences between the library angle and 
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(b) (c) 



FIG. 6: Data-MC comparison of fraction of energy in each of 
the 49 Csl crystals in an electron shower, (a) The fraction 
of energy in each of the 49 Csl crystals in electron showers 
for data, (b) KTeV03 data/MC ratio, (c) Current data/MC 
ratio. 



the desired angle are accounted for by shifting the 
incident position. The particle energy cutoff ap- 
plied in the GEANT shower library generation has 
been lowered from 600 keV to 50 keV for electrons; 
the photon cutoff of 50 keV is unchanged. Sixteen 
showers per bin have been generated. Energy de- 
posits are corrected for energy lost in the 12 /mi 
mylar wrapping around the Csl crystals. 

The current Monte Carlo produces a significantly 
better simulation of electromagnetic showers in the 
Csl calorimeter. Figure |3] shows the data-MC com- 
parison of the fraction of energy in each of the 49 
Csl crystals in electron showers relative to the to- 
tal reconstructed shower energies for electrons from 
Kl — > 7r ± e T ^ decays. The majority of the energy is 
deposited in the central crystal since the Moliere ra- 
dius of Csl is 3.8 cm. These plots are made for 16-32 
GeV electrons with incident angles of 20-30 mrad; 
the quality of agreement is similar for other ener- 
gies and angles. The data-MC agreement improves 
dramatically as seen in Fig. [SJ This improvement 
in the modeling of electromagnetic shower shapes 
leads to important reductions in the systematic un- 
certainties associcated with the reconstruction of 
photon showers from K — > 7r°7r° decays (see Sec. 

Unci) . 

2. Ionization Energy Loss. In KTeV03, we did not 
include the effect of ionization energy losses for 
charged particles in the simulation. In the current 
simulation, we include the ionization loss for each 
volume of material in the detector. The total loss 
up to the surface of the Csl is less than 5 MeV. 
This is a very small effect for K — > tt + tt~ decays 



but it is important for low-energy electrons used in 
the calibration of the Csl calorimeter and affects 
converted photons from K — > tt q tt° decays. 

3. Bremsstrahlung. In KTeV03, the MC included elec- 
tron Bremsstrahlung in materials upstream of the 
analyzing magnet only. In the current analysis, the 
Bremsstrahlung rate and photon angle in each vol- 
ume of material in the detector are included in the 
simulation. 

4. Delta Rays. In the KTeV03 simulation, delta rays 
produced in a drift chamber cell deposited all of 
their energy in that cell. The MC now has a more 
complete treatment in which delta rays may scat- 
ter into adjacent cells of the drift chamber. High 
momentum delta rays are traced through the de- 
tector like any other particle and low momentum 
delta rays are simulated using a library created 
with GEANT4[24}. This treatment of delta rays 
improves our simulation of the distribution of ex- 
tra in-time hits in the drift chambers. 

5. Pion Interactions. The probability for pions to in- 
teract hadronically with material in the spectrom- 
eter is 0.7%; hadronic interactions in the spectrom- 
eter were not simulated in KTeV03. These events 
are now simulated using a GEANT-based library 
which contains a list of secondary particles pro- 
duced by each hadronic interaction. An average 
of nine secondary particles are produced per inter- 
action; these secondary particles are read in from 
the shower library and traced through the rest of 
the detector like any other particle. Events with 
pion interactions typically trigger the photon veto 
system and so do not pass selection criteria in the 
analysis. 

6. Fringe Fields. The simulation of fringe fields from 
the analysis magnet has been refined. Fringe fields 
from the analysis magnet inside the vacuum tank 
and between all four drift chambers have been 
measured and are now simulated. The maximum 
strength of the fringe field in these regions is com- 
parable to the earth's magnetic field. The fringe 
field simulation improves the MC description of the 
azimuthal dependance of the it + tt~ invariant mass 
for K —> tt + tt~ data. 

7. Position Resolution. The position resolution of the 
drift chambers is dependent upon position within 
the cell as shown in Fig. [7J In KTeV03, the reso- 
lution was treated as flat across the cell; the posi- 
tion dependence of the resolution is now simulated. 
The position dependence of the resolution is also 
considered in the analysis of K — > tt + it~ data; this 
is described in Sec. IIII Bl 

The position resolution of Csl calorimeter clusters 
in the MC is slightly worse than in data. To better 
match the resolutions, we artificially improve the 
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FIG. 7: Dependence of drift chamber position resolution on 
position within the cell. D is the distance from the sense wire. 
Crosses represent the measured central values and uncertain- 
ties of the resolution in bins of D and the line represents a 
polynomial fit to the data. This fit is used in the simula- 
tion to parameterize the position dependence of the position 
resolution. 



range [— A s , +A S ] includes 68.3% of the area of a Gaus- 
sian with mean s and width <r s . 

B. Charged Reconstruction and Systematics 

The K — > 7r+7r~ analysis consists primarily of the re- 
construction of tracks in the spectrometer; the vertices 
and momenta of the tracks are used to calculate kine- 
matic quantities describing the decay. The Csl calorime- 
ter is used for particle identification by comparing the 
reconstructed cluster energy to the measured track mo- 
mentum. The analysis requirements provide clean identi- 
fication of well-reconstructed K — > 7r + 7r~ events with lit- 
tle background contamination. The cuts are sufficiently 
loose to reduce systematics from modeling of resolution 
tails. The K — > tt + tt~ reconstruction and event selection 
are described in the following sections; more details of 
the analysis are found in [201 ] . 

1. K — > tt + tv" Reconstruction 



resolution of the MC by 9%. This is done by mov- 
ing the reconstructed position toward the generated 
particle position. 



III. DATA ANALYSIS 

The analysis is designed to identify K — > tttt decays 
while removing poorly reconstructed events that are dif- 
ficult to simulate, and to reject background. For each 
decay mode, the same requirements are applied to de- 
cays in the vacuum and regenerator beams, so that most 
systematic uncertainties cancel in the single ratios used 

to measure \r]^ / p\ 2 and \rj no /p\ 2 - The following sections 

describe the analysis and the associated systematic un- 
certainties in Re(e'/e). 

A. Common Features 

Although many details of the charged and neutral de- 
cay mode analyses arc different, several features are com- 
mon to reduce systematic uncertainties. We select an 
identical 40-160 GeV/c kaon momentum range for both 
the charged and neutral decay modes. We also use the 
same z-vertex range of 110-158 m from the target for each 
decay mode. To simplify the treatment of background 
from kaons that scatter in the regenerator, the veto re- 
quirements for the charged and neutral mode analyses 
are made as similar as possible. 

When discussing systematic uncertainties, we typically 
estimate a potential shift s ± cr s , where s is the shift in 
Re(e'/e) and a s is the statistical uncertainty on s. We 
assign a symmetric systematic error, A s , such that the 



The spectrometer reconstruction begins by finding 
tracks separately in the x- and y- views. Track segments 
are found in the two drift chambers upstream of the mag- 
net and the two drift chambers downstream of the mag- 
net; these segments are then extrapolated to the center 
of the magnet. We require the extrapolated segments to 
match within 6 mm at the magnet mid-plane to form a 
combined track; they typically match to within 0.5 mm. 
Each particle momentum is determined from the track 
bend-angle in the magnet and a map of the magnetic 
field. 

The process of finding track segments depends on the 
alignment and calibration of the drift chambers. For the 
current analysis, we made new measurements of the drift 
chamber sizes and rotations. The survey of the wire posi- 
tions used a large coordinate measurement machine with 
a camera and magnifying lens mounted on the end of a 
movable arm. The measured drift chamber size is about 
0.02% larger than the nominal value found by scaling 
the 6.35 mm "cell" size. The relative non-orthogonality 
between DC1 and DC2 is limited to ±30 /irad. The un- 
certainty in Re(e'/e) associated with the drift chamber 
alignment and calibration is 0.20xl0 -4 . The momen- 
tum measurement uses the known kaon mass as a con- 
straint; the 0.022 MeV/c 2 uncertainty in the kaon mass 
corresponds to an uncertainty in Re(e'/e) of 0.10xl0~ 4 . 

If two x-tracks and two y-tracks are found, the recon- 
struction continues by extrapolating both sets of tracks 
upstream to define vertices projected on the x-z and y- 
z planes. The difference between these two projections, 
Az vtx , is used to define a vertex-^ 2 , 

xltx = {AZvtx/<JAzf , (5) 

where cta 2 is the resolution of Az v t x - This resolution 
depends on momentum and opening angle, and accounts 
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for multiple scattering effects. The two x-tracks and two 
y-tracks are assumed to originate from a common vertex 

ifx2*x<ioo. 

To determine the full particle trajectory, the x and 
y tracks are matched to each other based on their pro- 
jections to the Csl; the projected track positions must 
match Csl cluster positions to within 7 cm. 

An event is assigned to the regenerator beam if the 
regenerator x-position has the same sign as the in- 
coordinate of the kaon trajectory at the downstream face 
of the regenerator; otherwise, the event is assigned to the 
vacuum beam. 

In KTeV03 [20j , the track segments were reconstructed 
assuming that the position resolution of the drift cham- 
bers does not depend on the hit position within a cham- 
ber cell. To check this assumption, a special data sample 
was collected with the magnetic field turned off. Three 
chambers are used to reconstruct straight tracks and 
these tracks are compared to the hits reconstructed in 
the fourth chamber. The resulting position resolution 
(see Fig. [7]) shows a significant dependence on the dis- 
tance between the track and the sense wire. Tracks pass- 
ing close to a sense wire have worse resolution because of 
the time separation of drift electrons reaching the sense 
wire. The variation of the resolution for tracks close to 
the boundary of a drift cell can be partly explained by 
the electric field configuration. The drift-time depen- 
dent resolution is included in the MC and is used for the 
track segment reconstruction in the current analysis. The 
new resolution measurement improves reconstruction of 
^-vertex and track momentum. For example, the width 
of the 7r + 7r~ invariant mass is reduced by ~ 14%; this 
improvement is more significant for higher momentum 
kaons where it reaches ~ 25%. 

The kaon decay vertex position and the momenta of 
the two tracks forming the vertex are used to calculate 
the 7r + 7r~ invariant mass, their energy, and p\, the sum 
of their momenta transverse to the beam direction. 



2. K —¥ n tv Selection 

The K —> tt + tt~ event selection begins with the three- 
level trigger during data taking. Level 1 uses hits in 
the trigger hodoscopes and the drift chambers to select 
events consistent with two charged particles coming from 
the decay of a kaon that did not undergo large angle scat- 
tering in the defining collimator or regenerator prior to 
the decay. Level 2 uses custom hit counting electron- 
ics and a track finding system to select events with two 
tracks from a common vertex. The vertex requirement 
at trigger level is loose compared to the selection crite- 
ria in the offline analysis. The inefficiencies of the Level 
1 and Level 2 triggers are studied using Kl — > 7r ± e =F ^ 
decays from minimum-bias triggers. The uncertainty in 
Re(e'/e) associated with Level 1 and Level 2 event selec- 
tion is 0.2 x 1CP 4 . The Level 3 software filter reconstructs 
two charged tracks and makes loose cuts on reconstructed 



mass and particle indentification. To measure the Level 3 
inefficiency of the K — > tt + tt~ trigger, we perform the full 
offline analysis on "random accepts," a prescaled subset 
of the K — > 7r + 7r~ trigger that has no Level 3 require- 
ment. We find that the bias in Re(e'/e) from the Level 
3 trigger inefficiency is (0.30 ± 0.12) xlO -4 . We correct 
Re(e'/e) for this bias and assign an uncertainty in i?e(e'/ e) 
of0.12xl0~ 4 . 

The offline selection criteria for K — > 7r + 7r~ decays 
are tighter than those imposed by the trigger. The 
K — > 7r + 7r~ analysis requirements and any associated 
systematic uncertainties in Re{e'/e) are described in the 
following paragraphs. 

We make a number of cuts on energy deposits in the 
veto detectors. The most important veto requirement 
are the muon veto cuts, which suppress background from 
Kl —¥ 7r ± /i T j/ M decays, and the regenerator cuts, which 
reduce background from scattered kaons. Additional veto 
cuts are made for consistency with the K — > 7r°7r° anal- 
ysis. 

We also use the spectrometer and the calorimeter as 
"veto detectors." We reject events with any tracks other 
than those from the vertex. We require the ratio of re- 
constructed cluster energy to track momentum, E/p, to 
be less than 0.85 to identify the tracks as pions. We re- 
quire that the track momentum be greater than 8 GeV/c 
to ensure 100% efficiency for the muon veto detectors. 
These cuts suppress background from K e3 and de- 
cay modes. 

We remove events with 1.112 GeV/c 2 < m p7r < 
1.119 GeV/c 2 , where m p7r is the reconstructed invariant 
mass assuming the higher momentum particle is a proton. 
This removes background from A — > pir~ and A — > pir + 
decays where the proton is mis-identified as a pion. 

Figure [5] shows the 7r + 7r~ invariant mass dis- 
tributions for the two beams; we require 488 
MeV/c 2 < m 7r+7r - < 508 MeV/c 2 . Figure M shows the 
py distributions; we require p\ < 250 MeV 2 /c 2 . The pf, 
requirement is the only K — > tt + tt~ selection criterion 
that results in a systematic uncertainty in Re(e'/e). We 
vary the p| cut from 125 MeV 2 /c 2 to 1000 MeV 2 /c 2 and 
assign a systematic uncertainty in Re(e'/e) of O.lOxlO -4 
based on the change in Re(e'/e). 

To reduce our sensitivity to details of the Monte Carlo 
simulation, we require track trajectories to be clear of 
a number of physical apertures. We require that tracks 
point at least 2 mm into the Csl calorimeter away from 
the edges of the Collar Anti detector that surrounds the 
beamholes and at least 2.9 cm inside the outer edge of 
the Csl calorimeter. If the vertex position is upstream 
of the Mask Anti (MA, see Fig. [2]), we require that the 
track position at the MA be less than 4 cm in x and y 
from the nominal beam center. We cut away from wires 
at the edges of the drift chambers. To reduce the possi- 
bility of x and y track candidate mismatches, we require 
that the projections of the tracks at the Csl calorimeter 
be separated by 6 cm in x and 3 cm in y. We require that 
decays originate from within one of the beams by requir- 



9 



10 
10 ! 
10 : 
10' 
10 : 
10 ; 
10 

1 



• Data 
- MC 
MC Sig 



460 480 500 520 540 
MeV 2 /c 2 
Vac beam it + % mass 




460 480 500 520 540 
MeV/c 2 
Reg beam n + n mass 



FIG. 8: 7r + 7r~ invariant mass distribution for K — > tt + it~ 
candidate events in the vacuum (left) and regenerator (right) 
beams. The data distribution is shown as dots, the K — > 
7r + 7r~(7) signal MC (MC Sig) is shown as dotted histogram 
and the sum of signal and background MC is shown as a solid 
histogram. 
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FIG. 9: Pt distribution for K — > tv + ty~ candidate events in 
the vacuum (left) and regenerator (right) beams. The data 
distribution is shown as dots, the K —¥ ir + ir~( r y) signal MC 
(MC Sig) is shown as dotted histogram and the sum of signal 
and background MC is shown as a solid histogram. 



ing that the projection of the vertex (x,y) position along 
the kaon direction reconstructs inside a 75 cm 2 square at 
the z of the downstream edge of the regenerator. 

We require a minimum separation between the tracks 
in the x and y views at each drift chamber. This cut is 
defined in terms of the DC cell through which the track 
passes; we require that the tracks be separated by at least 
3 cells at each chamber. This track separation cut forms 
a limiting inner aperture and depends on the position of 
each wire within the drift chambers. The wire spacing is 
known with an uncertainty of 20 /im. There are varia- 
tions in the actual wire spacing, which are measured in 
data, but are not simulated in the Monte Carlo. To deter- 



mine the effect of these variations, we convolve the track 
illumination with the wire-cell size to determine the num- 
ber of events that migrate across the track separation cut 
in data but not in MC. We find that the bias in i?e(e'/e) 
is (-0.16 ± 0.12)xl0 -4 ; the corresponding uncertainty in 
Re(e'/e) is ±0.22xlO~ 4 

The effective regenerator edge, shown in Fig. G2 defines 
the upstream edge of acceptance for K —> it + tt~ decays 
in the regenerator beam. We find the effective regenera- 
tor edge by calculating the probability for two minimum 
ionizing pions to escape the last piece of scintillator with- 
out depositing enough energy to be vetoed. This calcu- 
lation depends on the measured average energy deposit 
of a muon passing through the regenerator Pb module, 
the fraction of energy coming from the last piece of scin- 
tillator due to the geometry of the phototube placement 
on the Pb module, the value of the trigger threshold, 
and the value of the offline cut on the energy deposit in 
the Pb module. We find the effective regenerator edge 
to be (1.65 ± 0.4) mm upstream of the physical edge 
in 1997 and (0.7 ± 0.4) mm upstream of the physical 
edge in 1999. The difference in effective edges is due to 
different offline cuts on the energy deposit in the Pb mod- 
ule. In 1997, the edge is defined by the trigger threshold. 
In 1999, a tight offline cut is applied. We evaluate the 
uncertainty in this measurement by varying the trigger 
threshold and the fraction of energy coming from the last 
piece of scintillator by ~15% each. The 0.4 mm uncer- 
tainty in the position of the effective regenerator edge 
leads to an uncertainty in Re(e'/e) of 0.20xl0~ 4 . 

We estimate the systematic uncertainty in Re{e'/e) as- 
sociated with the Monte Carlo simulation of drift cham- 
ber efficiencies by generating separate sets of MC in 
which scattering, DC efficiency maps, and accidental ac- 
tivity are turned off. We take 10% of the resulting vari- 
ation in i?e(e'/e), 0.15xl0 -4 , to be the systematic un- 
certainty associated with the simulation of drift chamber 
efficiencies. We vary the simulated drift chamber resolu- 
tions by 5% and, from the resulting variation in Re(e'/e), 
we assign a systematic error of 0.15xl0~ 4 . 

The systematic uncertainties in Re(e'/e) associated 
with the K — > 7r + 7r~ analysis are summarized in Ta- 
ble HI The total systematic uncertainty associated with 
the K — > 7r + 7r~ analysis is 0.81 xlO -4 ; this is reduced by 
~ 35% from KTeV03. 



C. Neutral Reconstruction and Systematics 

To reconstruct K 7r°7r° decays, we first identify four 
clusters of energy in the calorimeter and reconstruct the 
energies and positions of the photons associated with 
each cluster. A number of corrections are then made 
to the measured cluster energies based on our knowledge 
of the Csl calorimeter performance and the reconstruc- 
tion algorithm. We use the cluster positions and ener- 
gies along with the well-known pion mass to determine 
which pair of photons is associated with which neutral 
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Source 


Error on ReU'/e) (xlCT 


4 ) 




KTeV03 Result 


Current Result 






1997 1999 


Total 


LI and L2 Trigger 


0.20 


0.20 0.20 


0.20 


L3 Trigger 


0.54 


0.20 0.14 


0.12 


Alignment and Calibration 


0.28 


0.20 


0.20 


Momentum scale 


0.16 


0.10 


0.10 


Pt 


0.25 


0.10 


0.10 


DC efficiency modeling 


0.37 


0.15 


0.15 


DC resolution modeling 


0.15 


0.15 


0.15 


Background 


0.20 


0.20 


0.20 


Wire Spacing 


0.22 


0.22 


0.22 


Reg Edge 


0.20 


0.20 0.20 


0.20 


Acceptance 


0.79 


0.87 0.25 


0.41 


Upstream z 




0.33 0.48 


0.40 


Monte Carlo Statistics 


0.41 


0.28 0.28 


0.20 


Total 


1.26 


1.12 0.82 


0.81 



TABLE I: Summary of systematic uncertainties in Re(e'/e) 
from the K —¥ 7r + 7r~ analysis. For errors which are evaluated 
individually for each year, the individual errors are listed in 
columns and the total is the weighted average of the individ- 
ual errors. For those errors which are evaluated for the full 
dataset or taken to be the same for both years, only one num- 
ber is listed. The value of each systematic uncertainty from 
KTeV03 is provided for reference. 



pion from the kaon decay and to calculate the decay ver- 
tex, the center-of-energy, and the 7r°7r° invariant mass. 
The precision of the Csl calorimeter energy and position 
reconstruction is crucial to the K — > 7r°7r° analysis and 
has been improved significantly since KTeV03. Section 
IIII C II gives details of the Csl calorimeter reconstruction, 
Sec. IIII C 21 describes the reconstruction of K — > tt°tt 
decays, Sec. IIII C 31 describes the selection criteria for 
K —> 7r°7r° decays, and Sec. IIII C 41 describes the sys- 
tematic uncertainties associated with the Csl calorimeter 
energy reconstruction. 



1. Csl Calorimeter Energy and Position Reconstruction 

The first step in reconstructing clusters is to determine 
the energy deposited in each crystal of the Csl calorime- 
ter. We convert the digitized information to an energy us- 
ing constants for each channel that are determined from 
the electron calibration. An in-situ laser, which deliv- 
ers light at known intensities via quartz fibers to each 
Csl crystal, is used to calibrate the DPMTs and to mea- 
sure the less than 1% spill-to-spill drifts in each channel's 
gain. The "laser correction" removes these spill-to-spill 
changes and is applied before any clustering is performed. 

We define a "cluster" as a 7x7 array of small crystals 
or a 3x 3 array of large crystals. Clusters near the bound- 
ary between the small and large crystals (see Fig. [5]) may 
contain both sizes of crystals; in this case the cluster is 
defined as a 3 x 3 array of "large" crystals where the en- 
ergy deposit in four small crystals is summed to form a 



"large" crystal as needed. Each cluster is centered on a 
"seed crystal," containing the maximum energy deposit 
among the crystals in the cluster. An initial approxi- 
mation of the cluster energy is found by summing the 
energies of the crystals in the cluster. 

The x, y position of a cluster is reconstructed by calcu- 
lating the fraction of energy in neighboring columns and 
rows of crystals in the cluster. The x, y position algo- 
rithm uses a map that is based on assuming a uniform 
photon illumination across each crystal to convert these 
ratios to a position within the seed crystal. The position 
maps are made using isolated clusters from K — > tt q tt° 
data; no corrections to the position are applied based on 
incident particle angle. The final position is evaluated 
after all energy corrections are applied. 

The raw cluster energy must be corrected for a num- 
ber of geometric and detector effects. We apply "crystal- 
level" corrections that adjust the energy in each crys- 
tal that makes up the cluster and "cluster-level" correc- 
tions, which are multiplicative corrections to the total 
cluster energy. Many of the crystal-level corrections rely 
on "transverse energy maps" ; as a function of position 
within the seed crystal, these maps predict the distribu- 
tion of energy among the crystals within a cluster . They 
are made using isolated photon clusters from K — > 7r°7r° 
data. The crystal-level and cluster-level energy correc- 
tions are enumerated below. 

1. Partial Clusters. We correct for energy that is miss- 
ing from the cluster because of crystal energies that 
are below the readout threshold or because portions 
of the 3x3 or 7x7 cluster are located in the beam 
holes or outside the calorimeter. The energy in 
missing crystals is estimated using the transverse 
energy maps. The energy in crystals that were be- 
low threshold is estimated by a parameterization, 
which was determined from data, of the ratio of 
energy in a crystal to the readout threshold. The 
fraction of the readout threshold energy predicted 
to be present in a crystal decreases with distance 
from the seed crystal and increases logarithmically 
with cluster energy. 

2. Out of Cone. The "out-of-cone correction" is 
applied because an electromagnetic shower is not 
fully contained by the 7x7 small-crystal or 3x3 
large-crystal clusters. We determine the out-of- 
cone correction using the same GEANT simulation 
used to generate the Monte Carlo shower library 
(see Sec. Ill C|) . The correction is parameterized by 
a quadratic function of the reconstructed distance 
of the cluster position from the center of the seed 
crystal and a linear function of the reconstructed 
energy. The size of the correction varies by about 
1% across the face of a crystal and by about 0.2% 
per 100 GeV. There is no explicit dependence of the 
out-of-cone correction on incident angle; because 
the reconstructed positions are not corrected for 
incident particle angle, the angle effect is included 
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implicitly in our parameterization as a function of 
reconstructed position. The correction is generated 
separately for photons and electrons, and for small 
and large crystals. In KTeV03, the out-of-cone cor- 
rection was determined for small and large crystals 
using 8 GeV GEANT showers, but there was no 
adjustment for the energy, position, or type of the 
incident particle. 

3. Longitudinal Response. We correct the energy in 
each crystal for the ~5% non-uniformity of re- 
sponse along the length of each Csl crystal. The 
longitudinal response of each Csl crystal is mea- 
sured in ten 5-cm z bins using cosmic ray muons 
that pass vertically through the Csl calorimeter. 
These muons are detected by a cosmic ray ho- 
doscope consisting of three sets of 3 m-long, over- 
lapping plastic scintillation counters placed above 
and below the Csl calorimeter. Typically the crys- 
tal response increases with z as the shower nears 
the PMT. The measured Csl response is convolved 
with a GEANT prediction of the shower's longi- 
tudinal distribution to correct the energy in each 
crystal. The GEANT shower profiles are generated 
separately for photons and electrons. There are in- 
dividual profiles for each crystal position within the 
cluster; they are binned in local position relative to 
the center of the seed crystal and in the same six 
logarithmic cluster energy bins used in the Monte 
Carlo (see Sec. Ill Cp . The mean shower depth for 
photons and electrons varies logarithmically with 
energy. These crystal-by-crystal shower profiles are 
a significant improvement to the longitudinal uni- 
formity correction; in the KTeV03 analysis, the uni- 
formity correction was applied at cluster level based 
only on a predicted average longitudinal energy dis- 
tribution for a whole shower. 

4. Shared Energy. For clusters that overlap, we must 
partition the energy in the shared crystals. The 
"overlap correction" separates the energy deposited 
in two or more clusters that share crystals by using 
the transverse energy maps to predict how much 
energy each particle contributed to the shared crys- 
tals. 

The "neighbor correction" estimates the amount of 
underlying energy in each crystal that comes from 
nearby clusters that are less than 50 cm away but 
outside the 3x3 or 7x7 cluster boundary. The cor- 
rection uses a 13x13 map to predict the energy 
contribution from neighboring clusters. This map 
is similar to the transverse energy maps but does 
not depend on position within the Csl crystal and 
is generated using GEANT rather than data. 

We correct clusters near the beam holes for ex- 
tra energy that comes from nearby clusters that 
do not share crystals but which leak energy across 
the beam holes. This correction uses maps made 
using electrons from Kl — > 7r ± e =F z/ data. 



5. Detector Effects. We correct for a number of detec- 
tor effects including the observed transverse non- 
uniformity of energies across each crystal that re- 
mains after the out-of-cone correction, the non- 
linearity of each channel with energy which re- 
mains after the longitudinal uniformity correction, 
and global time variations in the Csl calorimeter 
response. These corrections are measured using 
E/p of electrons from Kl —> -k e^v decays, and 
are applied multiplicatively to the total cluster en- 
ergy. The "transverse non-uniformity correction" 
is made by dividing each cluster seed crystal into a 
5x5 grid and measuring the cluster energies of elec- 
trons in each of these position bins. A multiplica- 
tive correction is applied to the total cluster energy 
based on the cluster's reconstructed position within 
the seed crystal. The correction is normalized such 
that the average correction over each crystal (25 
bins) is 1.0. The "channel- by-channel linearity cor- 
rection" removes the residual energy non-linearity. 
It is measured separately for each Csl calorimeter 
channel in data and Monte Carlo. The "spill-by- 
spill correction" is applied to correct for time vari- 
ations in the response of the calorimeter as a whole; 
it is measured and applied separately for each spill. 
Each of these corrections has a maximum magni- 
tude of less than 1%. 



6. Photon- Electron Differences. For the K — > 7r°7r° 
analysis we apply a "photon correction" that is de- 
signed to remove any residual differences between 
photons and the electrons that are used to calibrate 
the calorimeter. The correction is based on pho- 
tons from K — > 7r°7r° and Kj, — > 7r 7r 7r° decays. 
It is measured separately for 1996, 1997, and 1999 
in nine regions of the calorimeter by fitting each 
event for the photon energies applying six (four) 
kinematic constraints for 7r°7r°7r (7r°7r°). The de- 
tails of the kinematic fit are described in [25j . This 
correction is most important for photons with ener- 
gies below 20 GeV; its magnitude is less than 0.2%. 
The photon correction is new for the current analy- 
sis; no correction of photon-electron differences was 
applied in KTeV03. 



The quality of the calibration and the Csl calorime- 
ter performance is evaluated by analyzing electrons from 
the Kl — > ^e^v calibration sample with all corrections 
applied. The electron calibration for 1996, 1997, and 
1999 is based on 1.5 billion total electrons. Figure [TU1 
shows the E/p distribution and the energy resolution as 
a function of momentum of these electrons after all cor- 
rections. The final energy resolution of the calorimeter 
is a E /E ~ 2%/VE® 0.4%, where E is in GeV. 
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FIG. 10: K e 3 electrons after all corrections, (a) E/p for 1.5 
x 10 9 electrons, (b) Energy resolution. The fine curve shows 
the momentum resolution function that has been subtracted 
from the E/p resolution to find the energy resolution. 



2. K 

K — > ir°ir° and Kl — > 7r 7r°7r° events are fully recon- 
structed using the positions and energies of the four 
or six photon clusters in the Csl calorimeter. The 
Kl —> 7r°7T 7r reconstruction is almost identical to the 
K — > 7r°7r° reconstruction, but for simplicity this discus- 
sion will be in terms of the ir°ir° reconstruction. Using 
cluster energies and positions, we are able to reconstruct 
the z vertex of the kaon decay, the (x,y) components of 
the center-of-energy of the kaon, the kaon energy, and 
the 7r°7r° invariant mass. 

We must first determine which pairs of photons are as- 
sociated in the K — > 7r°7r° decay. For four photons, there 
are three possible pairings. For each pairing we calculate 
di2, the distance in z between the w° decay vertex and 
Zcsi, the mean shower depth in the Csl crystals. Using 
the pion mass as a constraint, in the small angle approx- 
imation, we find the distance for each pair of photons to 
be 



\JE\E2 

"12 ~ 7-12, 



(6) 



where r\i is the transverse distance between the two pho- 
tons at the Csl calorimeter. 

For each pairing, we compare the calculated distance 
for each candidate pion. In most cases, only the correct 
pairing will give a consistent distance for both pions. The 
consistency of the measured distance is quantified using 
the pairing chi-squared variable: 
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In Eq. dij is the calculated distance for each pion, d avg 
is the weighted average of the distance dij for both pions, 
and Oij is the energy dependent vertex resolution for each 
pion. We choose the pairing that gives the minimum 
value of x 2 a- Using Monte Carlo events, we find that this 
procedure selects the wrong pairing for less than 0.01% 



of K — > 7T°7r decays in the final event sample. The z 
vertex of the kaon decay is taken to be TjQsI - d a „ s for 
the best pairing. 

Wc find the center-of-energy of the kaon decay at the 
Csl calorimeter plane by weighting the position of each 
photon with its energy. The x and y components of the 
center-of-energy are 



Vc 



(8) 



where the sums are over all four photons. The center- 
of-energy is the point at which the kaon would have in- 
tercepted the plane of the Csl calorimeter if it had not 
decayed, so we can calculate the (x,y) position of the de- 
cay vertex by assuming it lies on the line between the 
target and the center-of-energy. The x coordinate of the 
kaon decay vertex is used to determine whether the kaon 
came from the regenerator or the vacuum beam. 

The 7r°7r° invariant mass is calculated from the coor- 
dinates of the kaon decay vertex and the four photon 
positions and energies. The kaon energy is calculated as 
the sum of the four photon energies. 



3. 



The K — > 7r°7r° event selection begins with three lev- 
els of trigger requirements during data taking. The 
K — > 7T°7r Level 1 trigger requires that the total en- 
ergy in the Csl calorimeter be greater than 30 GeV. The 
inefficiency in this trigger is studied using K — > 7r + 7r~7r 
decays from the K — > 7r + 7r~ trigger; the inefficiency at 
a given energy, E tota i, is the ratio of events with energy 
greater than E tota \ for which the Level 1 trigger bit is not 
set to the total number of events with energy greater than 
Etotai- We find inefficiencies ranging from (0.5-1.6) x 10 
in the 40-160 GeV energy range used for the analysis. 
The impact of this inefficiency is slightly different in the 
vacuum and regenerator beams because of small differ- 
ences in the energy distributions. The resulting bias in 
Re(e'/e) is less than 0.02xl0 -4 , which we assign as a 
systematic error. 

The Level 2 trigger requirement is based on the "hard- 
ware cluster counter" (HCC)[26]. The inefficiency in this 
trigger is measured using Kl —> 7r 7r°7r° decays from a 
trigger that has no Level 2 requirement. We reconstruct 
the Kl —> 7r 7r 7r° decays without any requirement on the 
HCC; the Level 2 inefficiency is the ratio of the number 
of events that do not meet the HCC requirement to the 
total number of events found in the offline reconstruc- 
tion. The bias in Re(e'/e) produced by this inefficiency 
is determined using K — > ir°ir° MC. The inefficiency is 
simulated to within 10% by the Monte Carlo, so we take 
10% of the measured bias as the systematic uncertainty 
in Re(e'/e). The total uncertainty in Re(e'/e) associated 
with the Level 2 trigger is 0.19xl0~ 4 . 

The inefficiency of the Level 3 K — > 7r°7r° trigger is 
studied using "random accepts," a prescaled subset of 
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the K — > ir°ir° trigger that has no Level 3 requirement. 
We find no statistically signficant bias in Re(e'/e) and 
quote an uncertainty of 0.07xl0 -4 in Re(e'/e) based on 
the statistical precision of the bias measurement. 

The offline selection criteria for the K ~ > 7r°7r° sam- 
ple are designed to select events that are cleanly recon- 
structed, to suppress background, and to select kinematic 
and fiducial regions appropriate for the KTeV detector. 
We evaluate the systematic uncertainty associated with 
each of the following requirements by loosening or remov- 
ing the cut and evaluating the change in Re(e'/e). 

The energy of each Csl calorimeter cluster is required 
to be greater than 3 GeV because the clustering correc- 
tions and MC simulation are not reliable at very low en- 
ergies. The minimum distance between the reconstructed 
positions of Csl calorimeter clusters is required to be 
greater than 7.5 cm because it is difficult to separate 
the energy deposits in two very close clusters. Clusters 
very near the beam holes are not as well reconstructed 
because of energy leakage across the beam holes and mul- 
tiple overlapping or nearby clusters. In KTeV03, the in- 
ner Csl aperture was defined by the Collar Anti (CA) 
detector; we now remove events with clusters having a 
seed crystal in the first ring of crystals around the beam 
holes. We do not find any systematic variation of Re(e'/e) 
with these requirements. 

The variable describing the quality of the photon pair- 
ing, x^.o (Eq- El, is required to be less than 50. This is a 
rather loose cut since more than 99% of K — > tt°tt° events 
passing all other cuts have x 2 values below 10. The pri- 
mary purpose of this cut is to reduce background from 
Kl —> 7T 7r 7r° events in which two of the photons escape 
the detector; in this case it is likely that the missing 
photons come from different pions causing the remaining 
photons to be paired incorrectly. The systematic un- 
certainty in Re(e'/e) associated with this requirement is 
0.14xl0~ 4 . 

The "shape chi-squared" variable, x 2 , is a measure of 
how well the transverse energy distribution of each Csl 
calorimeter cluster matches the expected distribution for 
a photon. This variable, which is not a true chi-squared 
because of correlations that are not considered, is calcu- 
lated by comparing the transverse energy distribution of 
each cluster to the transverse energy maps described in 
Sec. IIII C II The maximum value of x 2 f° r each event 
is required to be less than 48. The purpose of this cut 
is to remove background from Kl — > 7r°7r 7r° events in 
which two or more photons overlap in the Csl calorimeter 
and are reconstructed as a single cluster. In these cases 
the transverse distribution of energy would tend to be 
different from that of a single photon cluster. The sys- 
tematic uncertainty in Re(e'/e) associated with the shape 
chi-squared requirement is 0.15xl0~ 4 . 

We make a number of cuts on the veto detectors to re- 
duce background. We also use the calorimeter, spectrom- 
eter, and trigger hodoscope as "veto detectors" by cut- 
ting on extra clusters, tracks, and hits. There is no sys- 
tematic uncertainty associated with these requirements. 
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FIG. 11: K — >• tv n RING distributions for data and signal 
MC in the vacuum (left) and regenerator (right) beams. The 
dashed line indicates our cut. 



In the K — » ir + ir~ analysis we use p^ to remove events 
in which the kaon scatters in the collimator or the re- 
generator. This variable is not available for K — > 7r°7r° 
decays since we do not measure the photon angles, so we 
use the "ring number" variable to reject scattered kaon 
decays. Ring number is calculated using the center-of- 
energy of the reconstructed clusters, and is defined as 

RING = 40000 x Max( Ax 2 coe , Ay 2 coe ), (9) 

where Ax coe and Ay coe are the distances from the center- 
of-energy to the center of the closest beam hole. A change 
of ARING= 1 corresponds to an incremental area of 1 
cm 2 centered on the beam hole. Events with ring number 
less than 81 cm 2 should be from kaons decaying inside 
one of the two beams. Figure [11] shows the ring number 
distributions for both beams for data and Monte Carlo. 
The ring number is required to be less than 110 cm 2 ; the 
systematic uncertainty in Re(e'/e) associated with this 
requirement is 0.27xl0^ 4 . 

The limiting apertures for K — > 7r°7r° events are the 
Csl calorimeter inner aperture at the beamholes, the Csl 
calorimeter outer aperture, the upstream edge in each 
beam, and an effective inner aperture resulting from 
the 7.5 cm photon separation requirement at the Csl 
calorimeter. The Csl calorimeter inner and outer aper- 
tures are defined by rejecting events in which a photon 
hits the innermost or outermost ring of Csl crystals. The 
upstream aperture in the vacuum beam is defined by the 
Mask Anti and the upstream aperture in the regenerator 
beam is defined by the lead module at the downstream 
edge of the regenerator. The systematic errors associated 
with the precision of these apertures are discussed in the 
KTeV03 paper[2(| and have not changed; the individual 
values are listed in Table [TTJ The total systematic uncer- 
tainty in Re(e'/e) associated with limiting apertures in 
the K -> 7r°7r° analysis is 0.48 xl0~ 4 . 

Figure [T2l shows the reconstructed kaon mass distribu- 
tions for both beams for data and Monte Carlo. The mass 
is required to be 490 MeV/c 2 < m^o < 505 MeV/c 2 . 
The sidebands of the m^o^o distribution are almost ex- 
clusively K~l — > 7r°7r 7r background, with a small con- 
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FIG. 12: K — > 7T 7r hIttOtt distributions for data and signal 
MC in the vacuum (left) and regenerator (right) beams. The 
dashed lines indicate our cuts. 



tribution from events in which the photons have been 
mispaired. The peaking background at the kaon mass is 
from decays of kaons which scattered with non-zero an- 
gle in the regenerator and the defining collimators. More 
details on the background are given in Sect. IIII D[ 



4- Energy Systematics 




I I l i i l I i I rW, | | | | | | | | | | | | | | | | | | | | | 

124 124.5 125 125.5 126 126.5 127 



z vertex (m) 

FIG. 13: Regenerator beam K — > 7r°7r° z vertex distribution 
near the regenerator for 1999 data and Monte Carlo, (a) Un- 
corrected data, (b) Data with energy scale correction applied. 



The reconstruction of K — > 7r°7r° decays depends en- 
tirely on the reconstruction of energies and positions of 
photon showers in the Csl calorimeter. Reconstructed 
quantities may depend upon the absolute energy scale 
or the energy linearity of the Csl calorimeter. We ap- 
ply corrections that match the energy scale between data 
and Monte Carlo, and we assign systematic uncertainties 
based on any disagreement in either absolute energy scale 
or energy linearity between data and Monte Carlo. The 
procedures for matching the energy scale and evaluating 
the energy systematics are described in this section. 

The energy scale of the Csl calorimeter is set by the 
electron calibration, but there is a small, residual differ- 
ence in energy scale between data and Monte Carlo for 
K -> 7T°7r events. This difference is removed by adjust- 
ing the energy scale in data such that the sharp edge 
in the z vertex distribution at the regenerator matches 
between data and Monte Carlo, as shown in Fig. [13] 
The correction is determined by sliding finely binned 
K — » ir°Tr° data and Monte Carlo z vertex distributions 
in the regenerator beam past each other and using the 
Kolmogorov-Smirnov (KS) test to determine how much 
the data must be adjusted to best match the MC. The 
correction is binned in kaon energy in the same 10 GeV 
energy bins that are used to extract our results (see Sec. 
lIVBp . The same correction is applied to each cluster in 
an event. 

The final energy scale adjustment is shown as a func- 
tion of kaon energy in Fig. QjJ The average size of the 
^-vertex shift is ~2.5 cm. This corresponds to an aver- 
age energy correction of ~0.04%, compared to ~0.1% in 
KTeV03. As a result of improvements to the simulation 
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FIG. 14: Change in the final energy scale adjustment rel- 
ative to KTeV03. The dashed line represents no data-MC 
mismatch. The y axis on the right side of the plot shows the 
data-MC z vertex shift in meters. 



and reconstruction of clusters, the required energy scale 
adjustment in the current analysis is smaller and less de- 
pendent on kaon energy than in the KTeV03 analysis. 

This final energy scale adjustment ensures that the en- 
ergy scale matches between data and MC at the regener- 
ator edge, but we must check whether the data and MC 
energy scales remain matched for the full length of the 
decay volume. Any non-linearity would result in different 
effective energy scales at different decay points because of 
the correlation between the z vertex and kaon energy dis- 
tributions. We check the energy scale at the downstream 
end of the decay region by studying the z-vertex distri- 
bution of 7T 7r° pairs produced by hadronic interactions 
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in the vacuum window and other downstream detector 
elements in data and MC. To verify that this type of 
production has a comparable energy scale to K —> ir°ir°, 
we also study the z- vertex distribution of hadronic 7r°7r° 
pairs produced in the regenerator. The z-vertex distri- 
bution of regenerator hadronic events is Gaussian while 
the distribution of downstream events is more compli- 
cated, as described below. The methods for making the 
data-MC comparison in each case are described in the 
following paragraphs. 

We compare the Gaussian z-vertex distributions of 
hadronically produced regenerator events between data 
and MC by sliding the distributions past each other and 
using the chi-squared test. The average data-MC differ- 
ence is plotted in Fig. [T5] we find no significant data-MC 
mismatch in this sample. 

For the downstream hadronic events, we consider in- 
teractions in four separate detector volumes: the vacuum 
window, the upstream drift chamber, and the two helium 
bags surrounding the drift chamber. The production of 
7r°7r° pairs in each of these volumes is simulated sepa- 
rately; a fit is used to determine the relative contribution 
of each material, and to find the difference between the 
data and MC z-vertex distributions. The fit is performed 
separately for the 1996, 1997, and 1999 data samples. 
Figure [16] shows the z-vertex distributions of downstream 
hadronic 7r°7r° pairs for 1999 data and MC, before and 
after the Monte Carlo data are shifted by the measured 
1.06 cm data-MC difference. The z shifts measured for 
each year are plotted in Fig. [15j 

To convert these shifts to an uncertainty in i?e(e'/e), 
we consider a linearly varying energy scale distortion such 
that no adjustment is made at the regenerator edge and 
the z shift at the vacuum window is that measured by the 
hadronic downstream sample. This distortion is shown 
by the shaded region in Fig. [T5J We rule out energy 
scale distortions that vary non-linearly as a function of 
z vertex because they introduce data-MC discrepancies 
in other distributions. The systematic error on Re(e'/e) 
due to uncertainties in the K — > 7r°7r° energy scale is 
0.65xl0" 4 . 

To evaluate the effect of energy non-linearities on the 
reconstruction, we study the way the reconstructed kaon 
mass, which does not depend on the absolute energy 
scale, varies with reconstructed kaon energy, kaon z ver- 
tex, minimum cluster separation, and incident photon 
angle. Data-MC comparisons for these distributions for 
the 1999 data sample are shown in Fig. [T7] To mea- 
sure any bias resulting from the nonlinearities that cause 
the small data-MC differences seen in these distribu- 
tions, we investigate adjustments to the cluster energies 
that improve the agreement between data and MC in 
the plot of reconstructed kaon mass vs kaon energy. We 
find that a 0.1%/100 GeV distortion produces the best 
data-MC agreement for the 1997 and 1999 datasets. Fig- 
ure [TS] shows the improvement in data-MC agreement 
with this distortion applied to 1999 data. The 1996 
dataset has slightly larger non-linearities; we find that 
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FIG. 15: Energy scale tests at the regenerator and vacuum 
window. The difference between the reconstructed z posi- 
tions for data and MC is plotted for K —¥ events, and 
for hadronically produced 7r°7r° pairs at the regenerator and 
the downstream detector elements. The solid point at the re- 
generator edge is the K — > sample; there is no difference 
between data and MC by construction. The open point at the 
regenerator edge is the average shift of the hadronic regener- 
ator samples for all three datasets. The points at the vacuum 
window are the shifts for the downstream hadronic events for 
each dataset separately. The shaded region shows the range 
of data-MC shifts covered by the total systematic uncertainty 
from the energy scale. For reference, the data-MC shift at the 
vacuum window from KTeV03 is also plotted. 
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FIG. 16: z-vertex distributions of n tt pairs produced 
hadronically in downstream detector elements for 1999 data 
and MC. (a) Data (dots) and nominal MC (histogram), (b) 
Data (dots) and MC that is shifted 1.06 cm downstream to 
match the data (histogram). 
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FIG. 17: Comparisons of the reconstructed kaon mass vs 
(a) z-vertex, (b) kaon energy, (c) minimum cluster separa- 
tion, and (d) photon angle for 1999 data (circles) and MC 
(stars). The values plotted are the difference between the re- 
constructed kaon mass for each bin and the PDG kaon mass. 
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FIG. 18: Effect of 0.1%/T00 GeV distortion on M K vs E K 
for 1999 data. The values plotted are the difference between 
the reconstructed kaon mass for each bin and the PDG kaon 
mass. 



a 0.3%/100 GeV distortion produces the best data-MC 
agreement for this dataset. The data-MC agreement in 
the reconstructed kaon mass as a function of kaon energy 
has been significantly improved compared to KTeV03, 
where a 0.7%/100 GeV distortion was required. 

To evaluate the systematic error associated with these 
non-linearities, we apply the distortions to the data and 
find that Re(e'/e) changes by less than 0.2 xl0~ 4 for all 
three datasets. Properly weighting the three datasets, we 
find that the systematic uncertainty on Re(e'/e) due to 
energy non-linearities is 0.15 xlO -4 . 

The systematic uncertainties in Re(e'/e) from the K — > 
7r°7r° analysis are summarized in Table |TTJ The K — > 



7r°7r° analysis contributes an uncertainty in Re(e'/e) of 
1.55xl0" 4 , which is reduced by -23% from KTeV03. 



Source 


Error on ReU'/e) (xl(T 4 ) 






KTeV03 Result 


Current Result 






1996 


1997 


1999 


Total 


LI Trigger 


0.10 


01 


0.01 


0.03 


09 


L2 Trigger 


n i ^ 


90 


0.12 


0.23 


1 Q 


L3 Trigger 


u.uo 


90 


0.04 


0.05 


07 


Ring Number 






0.27 




0.27 


Pairing: y 2 






0.14 




0.14 


Shape x 2 


0.20 




0.15 




0.15 


Energy Nonlinearity 


0.66 


0.10 


0.10 


0.20 


0.15 


Energy Scale 


1.27 


0.45 


0.82 


0.59 


0.65 


Position Reconstruction 


0.35 




0.35 




0.35 


Background 


1.07 


1.14 


1.06 


1.07 


Csl Inner Aperture 


0.42 




0.42 




0.42 


MA Aperture 


0.18 




0.18 




0.18 


Reg Edge 


0.04 




0.04 




0.04 


Csl Size 


0.15 




0.15 




0.15 


Acceptance 


0.39 


0.48 


0.48 


MC Statistics 


0.40 


0.75 


0.37 


0.41 


0.25 


Total 


2.01 


1.69 


1.63 


1.56 


1.55 



TABLE II: Summary of systematic uncertainties in Re(e'/e) 
from the K — > 7r°7r° analysis. For errors which are evaluated 
individually for each year, the individual errors are listed in 
columns and the total is the weighted average of the individ- 
ual errors. For those errors which are evaluated for the full 
dataset or taken to be the same for all years, only one num- 
ber is listed. The value of each systematic uncertainty from 
KTeV03 is provided for reference. 



D. Background and Systematics 

Background to the K — > irir signal modes is simulated 
using the Monte Carlo, normalized to data outside the 
signal region, and subtracted. There are two categories 
of background in this analysis: scattered K — > hit events 
and non-7T7r background. We use decays from coher- 
ently regenerated kaons only; any kaons that scatter 
with non-zero angle in the regenerator are treated as 
background. This regenerator scattering background and 
background from kaons that scatter in the defining col- 
limators have the same momentum and distributions 
for both K — » ir + ir~ and K — > 7r°7r° decays. This back- 
ground can be identified using the reconstructed trans- 
verse momentum of the decay products in the charged de- 
cay mode. Therefore, the scattering background is small 
in the charged mode, and we may use K — > n + ir~ decays 
to tune the simulation of scattering background on which 
we must rely in the neutral mode. 

Non-7T7r background is present because of misidentifi- 
cation of high branching-ratio decay modes. The back- 
ground to K — > ?T + 7r~ decays comes from Kl ix^e^v 
and Kl — > tt* decay modes. The background to 
K — s- 7r°7r° decays comes from Kl — > 7T 7r°7r° decays 
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and hadronic interactions in the regenerator. The back- 
ground estimation procedure and the associated system- 
atic uncertainties are described in detail in [20l |. 

There is only one significant change to the background 
estimation procedure since KTcV03. Hadronic produc- 
tion of K* and A resonances via Kl + N — > K s + X and 
n + N A + X are now included in the K — > 7r + 7r~ re- 
generator beam background analysis; these background 
sources were not considered in the KTeV03 analysis. The 
incident neutron spectrum is assumed to be the same as 
that of the A baryon, which is measured in data. For Kg 
decays, both K ± tt t and ir°Ks,Ks — > tt + tt~ modes are 
simulated. The K s — > n Ks background is normalized 
using the transverse momentum side band in the regen- 
erator beam. The K* s — > if ± 7r =F and A — > p^n^ decays 
are normalized using mass sidebands in the regenerator 
beam reconstructed assuming the vertex is located at the 
regenerator edge. These two modes are seperated using 
the momentum asymmetry distribution of the secondary 
particles. The hadronic K* and A background samples 
have negligible contributions to the signal region after all 
selection cuts, but including them improves the descrip- 
tion of mass and pt sidebands. 

The background levels in K — > 7r + 7r~ are illustrated in 
Fig. [5]and Fig. |H1 and the background to K — > tt°tt may 
be seen in Fig. Ql] and Fig. [12J Background contributes 
less than 0.1% of K — > 7r + 7r~ events and about 1% of 
K — > tt°tt° events. Tables Hill and llVl contain summaries 
of all the background fractions for each dataset. There 
are some variations in background levels among the years 
due to differences in trigger and veto requirements. The 
systematic uncertainty in Re(e'/e) due to background is 
0.20xl(T 4 from K -> tt+tt" and 1.07xlCT 4 from K -> 
7r°7r°. 



Source 


Vacuum Beam 
1997 1999 


Regenerator Beam 
1997 1999 


Regenerator Scattering 
Collimator Scattering 
Kl n ± e^ F v 
K L —5- ^fi* v 


0.009% 0.008% 
0.032% 0.032% 
0.034% 0.030% 


0.073% 0.075% 
0.009% 0.008% 
0.001% 0.001% 
0.001% 0.001% 


Total Background 


0.074% 0.070% 


0.083% 0.085% 



TABLE III: Summary of K — > tt + -k background levels. 



E. Data Summary 

The numbers of events collected in each beam are sum- 
marized in Table [V] After all event selection require- 
ments are applied and background is subtracted, we have 
a total of 25 million vacuum beam K — > ir + ir~ decays and 
6 million vacuum beam K — > 7r°7r° decays. 



IV. ACCEPTANCE AND FITTING 
A. Acceptance Correction and Systematics 

We use the Monte Carlo simulation to estimate the ac- 
ceptance of the detector in momentum and z- vertex bins 
in each beam. We evaluate the quality of this simulation 
by comparing z-vertex distributions in the vacuum beam 
between data and Monte Carlo. To account for small dif- 
ferences in the energy spectrum between data and Monte 
Carlo, we reweight the distributions, using the same 10 
GeV/c momentum bins used by the fitter (see Section 
HVBp . by adjusting the number of MC events in each bin 
so that the data and MC kaon momentum distributions 
agree. We fit the data-MC ratio of z-vertex distributions 
to a line, and call the slope of this line, s, the acceptance 
"z-slope." We use this z-slope to evaluate the systematic 
error on Re(e'/e). 

A z-slope affects the value of Re(e'/e) by producing 
a bias between the regenerator and vacuum beams be- 
cause of the different z vertex distributions in the two 
beams. A good approximation of the bias on Re(e'/e) is 
sAz/6 where Az is the difference of the mean z values 
for the vacuum and regenerator beam z vertex distribu- 
tions. The factor of 6 converts the bias on the vacuum- 
regenerator beam ratio to a bias on Re(e'/e). The values 
of Az are 5.6 m for the K — > tt + h~ sample and 7.2 m 
for the K — > 7r°7r° sample. We use the measured bias on 
Re(e'/e) and the statistical error on that measurement to 
assign a systematic uncertainty in Re(e'/e). 

The z-slopes for the full dataset are shown in Fig. [TIB 
We use the 25 million vacuum beam K —> tt + it~ de- 
cays to measure the z-slope in the charged decay mode. 
The uncertainty in Re(e'/e) associated with this z-slope is 
0.41 xl0~ 4 . We assign an additional uncertainty of 0.40 
xlO -4 based on a Re(e'/e) fit which excludes K — > ir + ir~ 
decays from the region upstream of the MA. This region 
is very sensitive to the value of the MA apperture cut 
(see Sec. IIIIBp . and, since it lies upstream of the re- 
generator edge, there are no regenerator beam decays to 
compensate for this dependence. Combining these two 
uncertainties, the systematic uncertainty in Re(e'/e) as- 
sociated with the acceptance correction is 0.57xl0 -4 . 

We also measure the data-MC z-slope in the high 
statistics Kl —> it e^v decay mode and find a slope 
that is similar in magnitude to the systematic uncertainty 
from K — » 7r + 7r~. We do not use the Kl — > it^e^v z- 
slope to set the systematic error because it is sensitive to 
different detector effects and has different particle types 
in the final state than K — > it + tt~ . 

We use 88 million Kl 7r 7r°7r° decays to measure 
the z-slope in the neutral decay mode. This mode has the 
same type of particles in the final state as K — ► ir°ir° , and 
it is more sensitive than 7r°7r° to potential problems in 
the reconstruction due to close clusters, energy leakage at 
the Csl calorimeter edges, and low photon energies. We 
assign an uncertainty on Re(e'/e) from the neutral mode 
acceptance of 0.48 x 10 -4 based on the Kl 7r 7r°7r° z- 
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Source 


Vacuum Beam 
1996 1997 1999 


Regenerator Beam 
1996 1997 1999 


Regenerator Scattering 
Collimator Scattering 

l\ L — y 7T 7T 7T 

Photon Mispairing 
Hadronic Production 


0.288% 0.260% 0.258% 
U.lUz/o U.lzz/o U.lzU/o 
0.444% 0.220% 0.301% 
0.007% 0.007% 0.008% 
0.002% 0.001% 


1.107% 1.092% 1.081% 
U.Uol/o U.Uyo/o U.Uyl/o 
0.015% 0.006% 0.012% 
0.007% 0.008% 0.007% 
0.007% 0.007% 0.007% 


Total Background 


0.835% 0.603% 0.678% 


1.209% 1.197% 1.190% 



TABLE IV: Summary of K — > ir n background levels. Note that photon mispairing is not subtracted from the data and is 
not included in the total background sum. 





Vacuum Beam Regenerator Beam 


K - 


■¥ 7T + 7r 


25107242 43674208 


K - 


* 7T 7T 


5968198 10180175 



TABLE V: Summary of event totals after all selection criteria 
and background subtraction. 



slope. We also measure the z-slope in K —> 7r°7r° decays 
and find that the results are consistent with those from 
Kl — > tt°tt°tt decays. There is no significant change in 
Re(e'/e) when the region upstream of the MA is excluded 
in the neutral mode, so no additional systematic uncer- 
tainty is required. 



B. Fitting and Systematics 

The value of Re(e'/e) and other kaon parameters Am, 
T5, e , and Im(e'/e) are determined using a fitting pro- 
gram. The fitting procedure is to minimize x 2 between 
background subtracted data and a prediction function. 
The prediction function uses the detector acceptance de- 
termined with the Monte Carlo simulation. The fits are 
performed in 10 GeV/c kaon momentum bins. There is 
no z binning to determine Re(e'/e), while a z-binned fit is 
performed to measure the other kaon parameters. Uncer- 
tainties from the fitting procedure are mainly related to 
regenerator properties and the dependence of the result 
on external parameters. 

Neglecting the contribution from K$ produced at the 
target (called target- Ks), the number of K — > irir events 
in the vacuum beam for a given p, z is 



V™(p,z)~T(p)\r 1 \ 2 ex V [-- 



(10) 



where t — [z 



^itt-k/p is the measured proper time 



relative to decays at the regenerator edge, rj = r/^ = 

e + e' (i] = 7/00 = e — 2e') for charged (neutral) decays and 
J~(p) is the kaon flux. The fitting program includes the 
contribution of ta rget-ifg by using a phenomenological 
model for K°/K° production at the target and propagat- 
ing the kaon states up to the decay volume. The model of 
target-iCg production is checked by floating the K° / K° 
flux ratio in the fit. The fitted fraction of target- Ks de- 



viates from the model by (2.5 ± 1.6)%. The associated 
systematic uncertainty in Re(e'/e) is ±0.12xl0 -4 . 
The number of events in the regenerator beam is 

N™(p,z)~T{p)T reg { V ) 
x [|p(p)| 2 exp(-^r) +M 2 exp(-^) 
+ 2|p(p)|M cos (Amt. + 4> p {p) - 4> v ) exp (-^ 

"(11) 

where p(p) is the momentum-dependent coherent regen- 
eration amplitude, <j> p (p) = arg(p), l/r ave = (1/t s + 
l/ r i)/2 and T reg (j>) is the relative kaon flux transmis- 
sion in the regenerator beam. The prediction function 
accounts for decays inside the regenerator by using the 
effective regenerator edge (Fig. |3Jd) as the start of the 
decay region. 

The parameters from Eqs. llOllTl are determined as dis- 
cussed below. The kaon mix, Tip), is a free parameter for 
each of the twelve 10 GeV/c momentum bins. Separate 
kaon fluxes are allowed for charged and neutral decays to 
account for slight differences in the data samples, so there 
are a total of 12 x 2 = 24 free fit parameters to describe 
the kaon flux. The flux ratio of decays in the vacuum and 
regenerator beams is, however, the same in both charged 
and neutral decay modes. The 1996 K — > tt°tt data has 
no corresponding K — > 7r + 7r~ data, so it is possible that 
there could be small differences in the flux ratio between 
the two years which do not cancel in the fit. We assign 
an uncertainty in Re(e'/e) of ±0.03xl0 -4 from this pos- 
sibility, following [i~8l |. 

The relative kaon flux attenuation in the regenera- 
tor beam, T reg , results from the shadow absorber and 
the regenerator itself. The attenuation is measured di- 
rectly from data by comparing the rate of Kl — > 7r + 7r~7r° 
decays in the vacuum (N^: ) and regenerator (iV+~°) 
beams. As mentioned earlier, a dedicated trigger was in- 
troduced in 1999 to improve the statistical precision of 
this measurement; the improved measurement is applica- 
ble for all years of data taking. The attenuation is found 
to be a linear function of p for a kaon momentum range of 
40 - 160 GeV/c as shown in Fig. HOI For the 61.5 GeV/c 
average momentum of Kl — > 7r + 7r~7r° decays, the trans- 
mission is T = (7.771 ± 0.004)%. The slope of the mo- 
mentum dependence is or = (—3.5 ± 0.2) x 10~ 5 c/GeV 
where the errors represent total experimental uncertain- 
ties. The uncertainty on the momentum dependence of 
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FIG. 19: (a) Comparison of the vacuum beam z distributions for data (dots) and MC (histogram), (b) The data-to-MC ratios 
are fit to a line, and the z-slopes (see text) are shown. All distributions are for the full data sample used in this analysis. 



the regenerator attenuation corresponds to a 0.08x10 
uncertainty in Re(e'/e). 

The regeneration amplitude p is related to the differ- 
ence between the forward kaon-nucleon scattering ampli- 
tudes for K°, /(0), and K°, /(0), [13 : 

/_ = S «. (12) 
P 



The KTeV regenerator is composed mostly of plastic scin- 
tillator with two thin lead plates in the last module (see 
Fig [3^). The dominant source for regeneration is forward 
scattering in carbon with small additional contributions 
from hydrogen and lead. The regeneration in hydrogen 
and lead is fixed in the fit while the parameters describing 
carbon regeneration are allowed to float. 

For an isoscalar target (C12) and high kaon momen- 
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FIG. 20: Ratio of Kl — > n + iv~n decay rates in the regen- 
erator to vacuum beam as a function of a kaon momentum. 
The uncertainty on the slope is statistical uncertainty of the 
measurement. 



turn, /L can be approximated by a single Regge trajec- 
tory [281 ] and the magnitude of /_ varies with kaon mo- 
mentum as a power law, |/_(p)| ~ p a . The analyticity of 
the forward scattering amplitude relates the magnitude 
of f-(p) and its phase, arg(/_). The magnitude of |/_| 
at p = 70 GeV/c and the power law a are the two free pa- 
rameters describing regeneration in the fit. We estimate 
the systematic error from the analyticity assumption by 
allowing <j> p to deviate 0.25° from analyticity; the associ- 
ated uncertainty in Re(e'/e) is 0.07xl0~ 4 . 

For scattering off complex nuclei, the effects of nuclear 
screening corrections are important and the regeneration 
amplitude cannot be described by a single power law. We 
vary the screening models in the fit and find that the as- 
sociated Re(e'/e) uncertainty is 0.24xl0~ 4 . More details 
on the screening corrections and the analyticity relation 
for the scattering amplitude are given in Appendix [A] 

The K~l lifetime, T£, is taken from [2{|; the uncertainty 
in Re(e'/e) due to the uncertainty in this measurement 
is 0.01 xl0~ 4 . The values of Am and t$ are fixed to our 
measurements (Eq.[2T]) for the Re(e'/e) fit and are floated 
in the z-binned fit. The uncertainty in Re(e'/e) due to 
the values of Am and ts used in the fit is 0.11 xlO -4 . 

The systematic uncertainties in Re(e'/e) associated 
with fitting are summarized in Table IVT1 The total sys- 
tematic uncertainty in Re(e'/ e) from fitting is 0.31 x 10 -4 . 



V. RESULTS 

A. Measurement of Re(e'/e) 

In the KTeV fit for Re(e'/e), the inputs are the ob- 
served number of K — » ir + ir~ and K — > 7r°7r° decays 
in each of twelve 10 GeV/c momentum bins. The kaon 
fluxes for K — > ir + ir~ and K — > 7r°7r° in each momentum 
bin, the regeneration parameters, and Re(e'/e) are free 



Source 


Error on Re(e'/e) 




(xlO" 4 ) 


Regenerator transmission 


0.08 


Target-ifs 


0.12 


Am and ts 


0.11 


Regenerator screening 


0.24 


4> p (analyticity) 


0.07 


1996 Ks/K L flux ratio 


0.03 


r L 


0.01 


Total 


0.31 



TABLE VI: Summary of systematic uncertainties in Re(e'/e) 
associated with fitting. 



parameters. CPT symmetry is assumed by setting the 
phases </>_| and 0oo equal to the superweak phase. The 



final KTeV result is: 
Re(e'/e) 



[19.2 ± l.l(stat) ± l.8(syst)} x 10~ 4 
[19.2 ± 2.1] x 10" 4 . (13) 



The fit x 2 is X 2 ' l v — 22.9/21. The systematic uncertain- 
ties in Re(e'/e) are summarized in Table IVlTI This result 
corresponds to a particle- antiparticle partial decay rate 
asymmetry of 



T(K° 



-T{K° 



T(K° 



T(K° 



.2 ±0.6) x 10~ 6 . 

(14) 



Source 


Error on 

K -> 7T + 


Re(e'/e 

7T~ K 


(xlO-*) 


— > 7T 7T 


Trigger 


0.23 




0.20 


Csl cluster reconstruction 






0.75 


Track reconstruction 


0.22 






Selection efficiency 


0.23 




0.34 


Apertures 


0.30 




0.48 


Acceptance 


0.57 




0.48 


Backgrounds 


0.20 




1.07 


MC statistics 


0.20 




0.25 


Total 


0.81 




1.55 


Fitting 


0.31 


Total 


1.78 



TABLE VII: Summary of systematic uncertainties in Re(e'/e). 
See Tables U and [jj] for more details on the errors from the 



K 



7i 7v analyses, respectively. 



We perform several checks of our result by dividing 
the data into subsets and checking the consistency of the 
Re(e'/e) result in the various subsets. To check for time 
dependence, we break the data into eleven run ranges 
with roughly equal statistics. There are five run ranges 
in 1997 and six in 1999. Since the 1996 K -> 7r°7T° data 
does not have any corresponding K — > ir + ir~ data, we 
combine it with the neutral mode data in the first 1997 
run range. As shown in Fig. 1211 we find consistent results 
in all of the run ranges. The combined result for 1996 
and 1997 data is Re(e'/e) = [20.0 ± 1.7(stat)} x 10~ 4 , 
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FIG. 21: Re(e/e) in subsets of the data sample. All points 
are statistically independent. The dashed line indicates the 
value of Re(e'/e) for the full data sample. The 97a run range 



includes the 1996 K 



data. 



which can be compared directly to the KTeV03 result of 
Re(e'/e) = [20.7 ± 1.5(stat)] x 10~ 4 20]. The decrease 
in statistical precision is due primarily to the removal of 
K — » 7r°7T events with a cluster near one of the beam 
holes (see Sect. IIII C 3[) . 

In 1999 we took data at high and low proton beam in- 
tensity so we are able to check for dependence of Re(e'/e) 
on beam intensity. About 43% of the 1999 data were col- 
lected at low intensity (defined as less than 1.25 x 10 11 
protons/s) while 57% were collected at high intensity 
(greater than 1.25 x 10 11 protons/s). The average proton 
rate in the low intensity sample is ~ 1 x 10 11 protons/s 
and the average rate at high intensity is ~ 1.6 x 10 11 
protons/s. The average rate for 1996 and 1997 data is 
~ 1.5 x 10 11 protons/s. Figure |2"2"1 shows the 1999 Re(e'/e) 
result for low and high intensities; there is no evidence 
for intensity dependence of the result. 

Figure |2"21 also shows the value of Re(e'/e) for the sub- 
sets of data with the regenerator in the left or right beam. 
We find no variation of Re(e'/e) with regenerator posi- 
tion, which shows that there is no significant intensity 
difference between the two neutral beams and that there 
is no significant left-right asymmetry in the detector. 

There are several crosschecks of the K — > 7r + 7r~ sample 
for which we do not divide the K — > tt°tt° sample. We 
divide the K — > n + ir~ sample based on the polarity of 
the analysis magnet and whether the tracks bend inward 
or outward in the magnet. In each of these cases, the 
K — > 7r 7r° sample is common to both data points and 
the errors are estimated by the difference in quadrature 
between the subset error and the nominal error. Figure 
|2"21 shows the Re(e'/e) results for each of these subsets; 
they all show good agreement. The fit results for tracks 
that bend in or out in the magnet are both slightly larger 
than the nominal result; in this case other fit parameters 
have changed in each fit to allow the higher values of 
Re{e'/e). 

We check for dependence on kaon momentum by break- 
ing the data into twelve 10 GeV/c momentum bins. In 
these fits, we fix the power-law dependence of the regen- 



FIG. 22: Re(e'/e) consistency with beam intensity, regenera- 
tor position, magnet polarity, and track bend. The intensity 
subsets are for the 1999 data only. Reg-left and reg-right refer 
to the position of the regenerator beam in the detector. These 
subsets are for the full data sample. Mag+ and Mag— are the 
magnet polarity and in/out are the bend of the two tracks in 
the magnet. In each of these subsets the K — > sample is 
common to both fits; the errors are estimated by taking the 
quadrature difference with the error for the full dataset. The 
dashed lines indicate the value of Re(e'/e) in the appropriate 
full data sample. 
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FIG. 23: (a) Re(e'/e) and (b) |/_(70 GeV/c)| in 10 GeV/c 
momentum bins. All points are statistically independent. The 
dashed lines indicate the values for the full data sample. 



eration amplitude to the value found in the nominal fit. 
The free parameters are Re(e'/e), |/-(70 GeV/c)|, and 
the charged and neutral kaon fluxes. Figure [23] shows 
the values of Re(e'/e) and |/-(70 GeV/c)| for these fits. 
We see no evidence for dependence of the Re(e'/e) result 
on kaon momentum. 
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Re(e7e) 



E731 
NA31 
NA48 

KTEV 



7.4 + 5.9 
23.0 ± 6.5 
14.7+ 2.2 

19.2+ 2.1 



uncertainties are evaluated for the fit, accounting for cor- 
relations among the parameters. CPT invariance is im- 
posed a posteriori, including the total errors of the pa- 
rameters with their correlations, to obtain precise mea- 
surements of Am and ts- 

This approach allows a self-consistent analysis of the 
data with and without CPT constraints. The results are 
crosschecked following the procedure developed in fl6L [Tsl . 
[20| , in which separate fits for Am and ts were performed 
with CPT invariance imposed a priori. 



10 20 30 (x1(T) 

New World Ave. «h 16.8+1.4 



1. Measurement ofrs, Am, (j> t , and Im(e' /e) with no CPT 
constraint 



FIG. 24: Final KTeV result for Re(e'/e) and comparisons 
with prev ious measurements from E73l|lq]. NA31[17j], and 
NA48pJ|- The new world average for Re(e'/e) is also shown. 



Our value of Re(e'Je) is also consistent with other ex- 
perimental results [TH [13, EH. The weighted average 
of the new KTeV result with previous measurements is 
Re(e'/e) = [16.8 ± 1.4] x 10~ 4 ; see Fig. [24] The consis- 
tency probability of these results is 13%. 



B. Measurements of Other Kaon Parameters 

The regenerator beam decay distribution is sensitive 
to the the kaon parameters ts, Am, (f> t , and Im(e'/e) 
(see Eq. Hi"]) . These parameters can be measured by 
fitting the decay vertex distribution in the regenerator 
beam. The analysis of the kaon parameters follows the 
general procedure developed in [ll| EE , but uses a 
new method to fit for all parameters simultaneously and 
apply the CPT constraints a posteriori. 

The "z-binned" fit uses both charged and neutral mode 
data. The data are binned in 2 meter wide z-bins in the 
regenerator beam from 124 m to 158 m. A single z-bin 
from 110 m to 158 m is used in the vacuum beam. As 
in the nominal i?e(e'/e) fit, there are 24 free parameters 
for the charged and neutral kaon fluxes in each of twelve 
momentum bins and two free parameters for the regen- 
eration amplitude and phase. Two additional parame- 
ters are used to fit the effective regenerator edge for the 
charged and neutral data to account for smearing effects 
near the edge. 

The five kaon parameters, ts, Am, <j) e , Re(e'/e), and 
Jm(e'/e) are free parameters of the z-binned fit. The fit 
thus provides the most general description of the data 
with no requirement of CPT invariance. All systematic 



The z-binned fit results are 



ts 



[89.589 ± 0.042 stat ± 0.056 syst 
[89.589 ±0.070] x 10~ 12 s 
[5279.7 ±12.8 st at±14.7 syst 
[5279.7 ± 19.5] x 10 6 H/s 
[43.86 ±0.40 sta t±0.49 syst ] c 
[43.86 ±0.63]° 
[21.10 ±1.31 st at±3.17 syst ] 
[21.10 ±3.43] x 10~ 4 



x 10" 



x 10 6 h/s 



x 10- 



Re(e'/e) 

Im(e'/e) = [-17.20 ± 9.04 stat ± 18.06 syst ] x HT 4 , 
= [-17.20 ± 20.20] x 10~ 4 . 

(15) 

The fit x 2 is X 2 /^ = 425.4/(432 - 33) The systematic 
uncertainties are summarized in Table IVIIII The total 
uncertainty and the correlations among the parameters 
are evaluated following the procedure described in Ap- 
pendix D of [2(| , and are given in Table IPX! The cor- 
relations among these results are shown in Fig. [55] and 
Fig. H3 

Uncertainties from the charged mode are smaller than 
those from the neutral mode for the measurements of 
Am, ts, and <j) e . The measurement of these parame- 
ters is effectively a statistical average of the charged and 
neutral mode values, so it is dominated by the statisti- 
cally larger charged mode. The measurements of Re(e'/e) 
and 7m(e'/e) depend on the difference between the two 
modes, so the uncertainties from the statistically smaller 
neutral dataset are more important. Similarly, the fit- 
ting uncertainties have a much larger impact on Am, rg, 
and <p<L than on Re(e'/e) and Im(e'/e) since the uncertain- 
ties in the regeneration properties enter directly for the 
former and cancel for the latter. 

There is a large correlation among ts, Am, and <p £ , 
and also between Re(e'/e) and Im(e' /e). The correlation 
between Am and e (p — 82.8%) is somewhat reduced 
compared to the pure statistical correlation (p — 97.3%), 
because the systematic uncertainty due to analyticity af- 
fects (/> e but not Am. 

The measurement of Im(e'/e) can be expressed in 
terms of the phase difference 



Ac/) w -37m(e'/e) 



[0.30±0.15 sta t±0.31 syst ] 
[0.30 ±0.35]°. 



(16) 
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It is consistent with zero as expected from CPT invari- 

ance in a decay amplitude. The individual values of </>_| 

and 0oo are 



b+- = [43.76 ±0.64]° 
6 00 = [44.06 ±0.68]°, 



(17) 



where the errors correspond to the total uncertainty and 
are calculated including the correlations reported in Ta- 
bleEO 

The superweak phase calculated using parameters from 
Eq.HUis 



<j) S w = [43.461 ± 0.069 sta t ± 0.070 syst ] 
= [43.461 ±0.098]°. 



The difference between d> f and 



>sw, 



[0.40±0.37 stat ±0.42 syst ] 
[0.40 ±0.56]°, 



(18) 



(19) 



is also consistent with zero as expected from CPT invari- 
ance in K°-K° mixing. 



2. Measurement of ts and Am with CPT constraint 

Large correlations of Am and t$ with e increase the 
experimental uncertainties on these parameters. Assum- 
ing 



4>t = 4>sw, 

Im(e'/e) = 0, 



(20) 



as required by CPT invariance, significantly reduces the 
errors. This effect is illustrated in Fig. [23 which shows 
Ay 2 contours of total uncertainty for Am, ts, and 4> e 
with and without the CPT constraint. 
The results are 



Am | cpt 

TS |cpt 

P 



= [5269.9 ± 12.3] x 10 6 ft/s 
= [89.623 ±0.047] x 10~ 12 
= -67.0%, 



(21) 



where the errors correspond to the total experimental 
uncertainty and p is the correlation coefficient between 
Am and ts ■ Compared to the determination without the 
CPT constraint, the uncertainty in Am is reduced by a 
factor of ~ 1.5. Using these values of Am and ts, we 
determine 



TABLE VIII: Systematic uncertainties for the global z-binned 
fit 
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n HQ 
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Track reconstruction 












maps 


a nnn 
U.UUU 


U.U 


n no 
U.UU 


U.U4 


A A Q 
U.45 


resolution 


0.001 


2.6 


0.08 


0.10 


1.20 


pt kick 


u.uuy 


U. i 


U.UU 


U. 14 


I. Id 


Z DC 


0.002 


0.1 


0.00 


0.28 


0.39 


Selection efficiency 












pt cut 


0.008 


3.6 


0.10 


0.16 


0.96 


accidental 


0.000 


0.1 


0.02 


0.05 


0.73 


scattering 


0.001 


0.3 


0.10 


0.15 


0.17 


Apertures 












Cell separation 


0.036 


10.0 


0.31 


0.42 


2.57 


Background 


0.001 


0.0 


0.01 


0.1 


0.6 


Acceptance 












Z slope 


0.007 


1.4 


0.04 


0.13 


3.05 


Trigger 


n nno 
U.UUz 


u . y 


n no 
U.Uz 


U.Uo 


1 71 
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Csl Reconstruction 












Energy linearity 


U.UUo 
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O Qfl 
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O A Q 


Energy scale 


U.UUo 


U.o 


U.U1 
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Selection Efficiency 












Ring 


0.002 


0.3 


0.01 


0.18 


2.19 


Pairing v 2 


0.012 


2.2 


0.07 


0.02 


2.19 


Shape \ 2 


0.0 


0.2 


0.02 


0.06 


0.90 


Apertures 












Csl size 


0.006 


0.2 


0.04 


0.64 


8.35 


MA 


0. 


0.1 


0.00 


0.27 


0.21 


CA 


0. 


0.2 


0.01 


0.47 


0.32 


Background 


0.008 


0.3 


0.04 


0.43 


6.69 


Acceptance 


0.002 


0.1 


0.01 


0.13 


2.81 


Fitting 












Attenuation Norm 


0.003 


0.3 


0.01 


0.01 


0.01 


Attenuation Slope 


0.003 


2.1 


0.05 


0.05 


0.00 


Target K s 


0.026 


4.7 


0.11 


0.00 


0.00 


Screening 


0.018 


5.6 


0.02 


0.57 


1.35 


Analytisity 


0.0 


0.0 


0.25 


0.0 


0.0 


MC statistics 


0.016 


4.9 


0.15 


0.36 


2.78 


Total Syst 


0.056 


14.7 


0.49 


3.17 


18.06 


Stat Error 


0.042 


12.8 


0.40 


1.31 


9.04 


Total Error 


0.070 


19.5 


0.63 


3.43 


20.20 



TABLE IX: Correlation coefficients for the z-binned fit. 



fowlcpt = [43.419 ±0.058]°. (22) 
3. Kaon Parameter Crosschecks 



ts Am (f> t Re(e'/e) 
Am -85.8% 
<j) t -76.5% 82.8% 

Re{e'/e) 6.8% -6.6% 2.5% 
Im(e'/e) -1.0% 2.6% -4.1% -46.5% 



We compare the new procedure to determine Am \ cpt 
and ts \ C pt to the one used in [20| in several steps. First 
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FIG. 25: A\ 2 = 1 contours of total uncertainty for (a) Am- 
ts, (b) <j} e -Am and (c) rs-cpe- Larger ellipses correspond to 
the 2-binned fit without CPT invariance assumption. Dashed 
lines correspond to (f> t = 4>sw CPT constraint. Smaller el- 
lipses are obtained after applying this constraint. 
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FIG. 26: Ax 2 = 1 contour for Re(e'/e) vs Im(e'/e) as mea- 
sured by KTeV compared to the measurement of -kit phase 
shifts [2j and the CPT invariance expectation. 



the CPT constraint is applied using statistical uncer- 
tainties only; these results agree exactly with the fits 
in which the CPT constraints are included as in [2fj| . 
The value of Am obtained using statistical uncertain- 
ties only, Am | s tat. cpt = 5266. 5xl0 6 is somewhat 
lower than for the determination using full uncertainties, 
Eq. [21] This difference can be traced to the above men- 
tioned reduction in the correlation between Am and </> e 
when systematic uncertainties are included. 

Next, Am | cp t and ts \ C pt are determined in the neu- 
tral and charged modes separately, as was done in 
KTeV03[13. The resulting values are Am 00 = (5 2 5 7.6 ± 
8.3) x 10 6 ^/s, Am+~ = (5269.0 ± 4.2) x 10 6 ?i/s, r s 00 = 



(89.667±0.039)xl0 

-12 



-12 



s, and Ts 



(89.620±0.020)x 

10~ 12 s, where the superscript H — (00) stands for the 
charged (neutral) mode and the errors represent statisti- 



cal uncertainties only. The measurements agree to within 
1.3cr st at- Finally, the total uncertainties from Eq. 1211 are 
compared to an evaluation in which the CPT constraints 
are embedded in the fit. They agree to within ~ 10%, 
which is consistent with small changes in correlations 
among the fit parameters. 



4- Determination of K°-K° Mass Difference 

KTeV measurements of the kaon system parameters 
can be u sed to determine the mass difference between K° 
and K°, which is zero in the absence of CPT violation. 
This test uses the Bell-Steinberger relation [3(|, which 
connects the CP and CPT violation in the mass matrix 
to the CP and CPT violation in the decay. Following the 
notation used in [3l|, the Bell-Steinberger relation can 
be written as 



r s + r L 



+ i tan 4>sw 



Re(e) 



Hm(S) 



r.« 



J2MfW s (f), (23) 
/ 



where the sum runs over all final states / and A^^if) = 
A(K L ,s — > /)■ The parameter S is related to the K°-K° 
mass and decay width difference: 



i (m K o — m^o) + \ (r^o — T^-o) 



i(f)sw 



■ cos <psw e 

1 S - 1 L 

For neutral kaons, only a few decay modes contribute 
significantly to the sum in Eq. [231 The largest contribu- 
tion comes from the Kl,s — > 7r + 7r~ and Kl,s — > 7r°7r° 
decay modes for which we can use KTeV measurements 
only; the other required measurements are the results 
from Eq. [151 and the Kl — > 7r + 7r~ and Kl — > 7r°7r° 
branching fraction measurements from KTeV 32]. The 
only external input needed is the value of tl ■ We use the 
PDG average, which is based mainly on measurements 

fromKLOE lHHl]. 

For the hadronic modes, we define 



1 



A L (i)A s (i) = Vi B(K s ^i), 



- 7r 7r°,7r+7r"(7),37r u ,7r u 7r + 7r"(7) 



(25) 



For the Kl — > t: + tt and Kl — > 7r°7r° decay modes we 
find 

cv^- = [(1124 ± 13) + i(1077± 13)] x 10~ 6 , (26) 
a^o = [(481 ± 7) -M(465 ± 7)] x 10~ 6 . (27) 



For the Kl —> 7r + 7r it and Kl 
use the PDG values 



37T° decay modes we 



a„+„- n o = [(0±2) + i(0±2)] x 10 
la^o^o^ol < 7 x 10~ 6 at 95% C.L.. 



-6 



(28) 
(29) 
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FIG. 27: Allowed region at 95% C.L. for AM, AI\ 



For the semileptonic decay modes we use the definition 
and values from the PDG 

a nlv = — y~) A L {irlv)A* s (irlv) + 
1 s , 

irlv 

+2i^B(K L -> Ttlv)Im{8) (30) 

and 

a„i„ = [(-2±5)+*(l±5)] x 10~ 6 . (31) 

Combining the experimental data, we find from Eq. [23] 

Im(S) = (-1.5 ± 1.6) x 10~ 5 . (32) 

The uncertainty in Im(S) is 1.9 x 10~ 5 in the PDG av- 
erage; the reduction of this error in Eq. Incomes mainly 
from the more precise measurement of </> £ — (f>sw from 
KTeV. 

Combining the new determination of Im{§) with the 
determination of Re(6) and Im(8) from the charge asym- 
metry in semileptonic decays (35j, 

Re(S) = (30 ±23) x 10~ 5 , 

Im(8) = (-660 ±650) x 10~ 5 , 

p = -21%, (33) 

we find 

Re{6) = (25.1 ±22.5) x 10~ 5 , 

Im{8) = (-1.5 ±1.6) X 10" 5 (34) 

with negligible correlation between the real and imagi- 
nary parts. 

Using the value of 6 from Eq. [MJ we derive the allowed 
region for AM = m R o — mg-o and Ar = T^-o — T^-o as 



shown in Fig. [27j Assuming no CPT violation is present 
in the decay amplitudes, i.e. Ar = 0, we obtain the limit 
for the mass difference, 

\M K o - M R0 \ < 4.8 x 10~ 19 GeV/c 2 at 95% C.L. (35) 



VI. CONCLUSIONS 

Using the full data sample of the KTeV experiment, 
we have made improved measurements of the direct CP 
violation parameter, Re(e'/e), and other parameters of 
the neutral kaon system. All of these results supersede 
previous KTeV results. 

Assuming CPT invariance, we measure the direct CP 
violation parameter 

Re(e'/e) = [19.2 ± l.l(stat) ± l.8{syst)} x 10" 4 (36) 
= [19.2 ±2.1] x 10~ 4 . 

Also under the assumption of CPT invariance, we report 
new measurements of the Kl, — K$ mass difference and 
the Ks lifetime: 

Am = [5269.9 ± 12.3] x 10 6 ft/s , , 

t s = [89.623 ±0.047] x 10" 12 s. ^ ' 

To test CPT symmetry, we measure the phase differ- 
ence between the K — > ir + ir~ and K — > 7r°7r° decays to 
be 

A0 = -3Im(e'/e) ( 
= [0.30 ±0.35]°, [M) 

and the phase difference relative to the superweak phase 
to be 

cj) e -<t> sw = [0.40 ±0.56]°. (39) 

These phase results are consistent with CPT invariance 
in both the decay amplitudes and K° — K° mixing. As- 
suming no CPT violation in the decay amplitudes, we set 
a limit on the K° — K° mass difference: 

\M K o - M R0 \ < 4.8 x 10" 19 GeV/c 2 at 95% C.L. (40) 

After decades of experimental effort, direct CP vio- 
lation in the neutral kaon system has now been mea- 
sured with an uncertainty of about 10%. Considerable 
improvement in theoretical calculations of Re(e'/e) will 
be required to take advantage of this experimental preci- 
sion. There is some optimism [13 ]. however, that future 
calculations using lattice gauge theory may approach a 
10% uncertainty, making the precise measurement of e'/e 
an equally precise test of the Standard Model. 



26 



Acknowledgments 

We gratefully acknowledge the support and effort of 
the Fermilab staff and the technical staffs of the par- 
ticipating institutions for their vital contributions. This 
work was supported in part by the U.S. Department of 
Energy, The National Science Foundation, The Ministry 
of Education and Science of Japan, Fundacao de Amparo 
a Pesquisa do Estado de Sao Paulo-FAPESP, Conselho 
Nacional de Desenvolvimento Cientifico e Tecnologico- 
CNPq and CAPES-Ministerio Educacao. 



Appendix A: Screening Corrections to Kaon 
Regeneration and Analyticity 

The measurements of (f> e and Am rely on an accurate 
determination of the regeneration phase. The regenera- 
tion phase is calculated from the momentum dependence 
of the regeneration amplitude using an analyticity rela- 
tion. Regeneration in carbon, the principal material of 
the KTeV regenerator [42j], is discussed here. We show 
how the momentum dependence of the regeneration am- 
plitude, which is calculated based on Regge theory with 
nuclear screening corrections, is checked using K — > 7T7t 
data. Based on this study, we derive the systematic un- 
certainties of the kaon parameters which are correlated 
with the determination of the regeneration phase. 

In Regge theory, for an isoscalar target and momenta 
above 10 GeV/c, the momentum dependence for the mag- 
nitude of the forward kaon-nucleon scattering amplitude 
difference for K° and K° (Eq. [T"2"|) is given by a single 
power law [28j : 



\f-(p)\ = \.f P (70 GeV/c)| 



P 



70 GeV/c 



(Al) 



Analyticity of the forward scattering amplitudes leads to 
a kaon momentum independent phase of /£" given by 



arg(/^; 



(A2) 



Therefore, a and |/£"(70 GeV/c)| determine f£(p) fully; 
they enter as free parameters in fits to KTeV data. 

Kaon-nucleon interactions in carbon are screened due 
to rescattering processes. The effect of screening modifies 
the momentum dependence of |/_(p)| as well as its phase: 



f-(p) = f P (p)-S sc (p)e 



i(f> sc (p) 



(A3) 



Here S sc (p) and 4> sc (p) stand for screening corrections, 
which are evaluated using Glauber theory formalism [36], 
\3l\ for diffractive scattering, and using various mod- 
els [HI for inelastic scattering. Elastic screening cor- 
rections are calculated with small theoretical uncertain- 
ties. Inelastic corrections, however, have a large spread 
of predictions as shown in Fig. 1281 Since a momentum 
independent change of |/_| does not modify arg(/_), cal- 
culations of the screening corrections fix 6 sc (p) = 1 for 
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FIG. 28: Kaon momentum dependence of screening correc- 
tions to (a:) the magnitude and (b:) the phase of f-(p), for 
various models from [3^]. The magnitude of each correction 
is defined to be 1 at 70 GeV/c. 
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FIG. 29: Kaon momentum dependence of (a:) ratio of 
from the fits using elastic (solid line) and inelastic factorized 
(dashed line) screening corrections to |/£"(p)| from the fit us- 
ing no screening corrections; (b:) arg(/_(p)) from the fits 
with elastic (solid line), inelastic (dashed line) screening cor- 
rections and with no screening corrections (dotted line). 



p = 70 GeV/c. Screening corrections calculated using the 
inelastic factorized model [38| are used for central values 
of all KTeV results in this paper. 

As shown in Fig. I28l screening corrections modify the 
value of effective power law a from Eq. IAH Since a 
is a free parameter of fits to the data, it is adjusted 
such that the net effect of the screening corrections is 
reduced. Figure l29h shows ratios of fitted [43| |/-(p)| 
using elastic and inelastic factorized screening models to 
fitted \ f£(p)\ using no screening corrections. Figure l2"5b 
shows the phase, arg(/_(p)), obtained in these fits 44j . 
The screening correction thus results in up to a 0.5% 
adjustment of the regeneration amplitude and up to a 
2° change in the regeneration phase, for the momentum 



27 



TABLE X: Power law coefficient a and \ 7 ' l v f° r fits using var- 
ious screening models. The uncertainties given in parenthesis 
are statistical errors. 





a X /v 


No screening 

Elastic screening 

Inelastic factorized screening 

Inelastic symmetric screening 


-0.5813(5) 471/399 
-0.5715(5) 431/399 
-0.5376(5) 425/399 
-0.5803(5) 438/399 



range of the analysis. The values of a obtained in the 
fits with different screening models are given in Table IXl 
Comparing y 2 /v for these fits, one can see that the data 
rule out calculations without screening corrections and 
disfavor the inelastic symmetric screening calculation. 

Comparison of Fig. [55] and Fig. suggests significant 
sensitivity to the regeneration parameters and that the 
screening corrections may be determined directly from a 
fit to the K —> hit data. In this fit, the magnitude of the 
regeneration amplitude, |/_(p)|, is floated as a free pa- 
rameter in each of the twelve 10 GeV/c momentum bins, 
and is approximated to be constant within each bin. The 
regeneration phase, arg(/_(p)), is calculated dynamically 
using the Derivative Analyticity Relation (DAR) (39| : 



arg(/-(p)) = -7T - W|^r 



)ln|/-(p)|. (A4) 



In this calculation, arg(/_(p)) is estimated at an aver- 
age of bin centers for the neighboring momentum bins, 
assuming that the bin-to-bin dependence of |/_(p)| is de- 
scribed by a power law. For momenta away from the bin 
centers, the phase is interpolated linearly. For example, 
the local power law at p = 50 GeV/c is calculated as 



In 



a(50 GeV/c) 



/_(55 GeV/c) 



/_(45 GeV/ C ; 



In 



55 



(A5) 



The phase arg(/_(50 GeV)) is then given by Eq. IA2I 
This fit is referred to as the "DAR fit". Compared to 
the standard z-binned fit, the DAR fit has an additional 
12—1 free parameters. 

Figure l3"0b compares ratios of |/_(p)| obtained from 
the DAR fit and from the fit using the inelastic factor- 
ized screening correction to | /_ (p) | obtained from the fit 
using no screening corrections. There is a good agree- 
ment between the inelastic screening calculation and the 
DAR fit. The quality of the DAR fit {^jv = 412/388) 
is similar to that of the fit using inelastic screening cor- 
rections. 

The data precision makes it possible to go one step 
further and float both |/_(p)| and arg(/_(p)) for each 
momentum bin, this fit is called the "phase fit" . An over- 
all phase can not be measured this way, since it is 100% 
correlated with the value of <f> e . Therefore the phase is 
fixed at p = 75 GeV/c to the value obtained in the DAR 
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FIG. 30: Kaon momentum dependence of (a:) ratio of |/_(p)| 
from the fit using inelastic factorized screening corrections 
(dashed line) and from the DAR fit (dots with error bars) 
to \f^(p)\ from the fit using no screening correction; (b:) 
arg(/_(p)) for the DAR fit (solid line), the phase fit (dots 
with error bars) and the fit using inelastic factorized screen- 
ing model (dashed line). The phase fit is fixed to the DAR 
value for 70 GeV/c < p < 80 GeV/c momentum bin as indi- 
cated by the arrow. 



fit. The phase fit thus does not check an absolute pre- 
diction of the DAR fit but only the predicted momentum 
dependence. The quality of the phase fit is not improved 
compared to the DAR fit: for extra 11 degrees of free- 
dom x 2 is reduced by 7 units only. Figure l3"0b compares 
arg(/_(p)) calculated in the DAR fit and in the fit us- 
ing the factorized inelastic screening correction as well 
as directly fitted in the phase fit. These three different 
determinations are consistent with each other. 

Note that the large correlation between arg(/_) ex- 
tracted at different momenta in the phase fit complicates 
a quantitative comparison between the predictions, so 
an additional test has been performed. The variation 
of arg(/_) versus momentum, predicted by the DAR fit, 
shows an approximately linear dependence. The total 
variation, for the 40 — 160 GeV/c momentum range, is 
about 3°. As a crosscheck, linear dependence of arg(/_) 
is assumed: 



arg(/-(p)) = arg(/ p ) + a R ■ 



p- 70 GeV/c 
120 GeV/c ' 



(A6) 



where p is in GeV/c and an is an additional free fit pa- 
rameter. A fit with |/_(p)| floated in twelve momentum 
bins using no screening corrections leads to significantly 
non-zero an: aji = (3.1 ± 0.6)°. A similar fit, which 
includes the DAR phase correction, gives an consistent 
with zero: a R = (0.2 ± 0.6)°. 

In KTeV03 [20], the systematic uncertainties due to 
the screening corrections were evaluated by comparing 
the factorized and symmetric inelastic screening models 
which led to a large 0.75° error for <p e . In this analysis, 
the systematic uncertainties are determined directly from 
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the data and estimated by comparing the fit based on the 
inelastic factorized screening calculation with the DAR 
fit. For (f) e , this evaluation leads to a small uncertainty 
of 0.02°, see Table ["VTITl The uncertainty is larger for 
Am (5.6 x 10 -12 h/s), since Am is more sensitive to the 
variation of arg(/_) as a function of the kaon momentum. 

The determination of the regeneration phase relies on 
the validity of the analyticity assumption, Eq. IA2I (or 
equivalently of Eq. IA4|) . Various sources of deviation 
from Eq. IA2I were studied in • The net effect of sub- 
leading Regge trajectories, uncertainties in /_ at low mo- 
mentum, and electromagnetic regeneration in lead were 
estimated to be below 0.25°, which is taken as a system- 
atic uncertainty for the regeneration phase. Note that 



an exchange of an Odderon can affect the regeneration 
amplitude at high momentum, but is not taken into ac- 
count in this evaluation. The Odderon contribution can 
be limited using pp and pp data; the best fit to these data 
gives a regeneration phase shift of —(0.2 ± 0.6)° [40]. 

A cross check of the analyticity assumption is per- 
formed by measuring the asymmetry between K — > 
■K~e + D e and K — > n + e~i> e decays downstream of the 
regenerator which is sensitive to the regeneration phase. 
Using approximately 125 million decays, the difference 
between the measured arg(/_) and the prediction based 
on analyticity is [-0.70 ± 0.88(stat) ± 0M(syst)}° [H|, 
which is consistent with the analyticity assumption. 
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Total regeneration amplitude of the KTeV regenerator 
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Note that small variation of the phase for the fit using 
no screening corrections is caused by the regeneration in 
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